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Abstract 

We derive a set of integral equations of the TBA type for the generahzed cusp anomalous 
dimension, or the quark antiquark potential on the three sphere, as a function of the angles. 
We do this by considering a family of local operators on a Wilson loop with charge L. In 
the large L limit the problem can be solved in terms of a certain boundary reflection matrix. 
We determine this reflection matrix by using the symmetries and the boundary crossing 
equation. The cusp is introduced through a relative rotation between the two boundaries. 
Then the TBA trick of exchanging space and time leads to an exact equation for all values 
of L. The L = case corresponds to the cusped Wilson loop with no operators inserted. We 
then derive a slightly simplified integral equation which describes the small angle limit. We 
solve this equation up to three loops in perturbation theory and match the results that were 
obtained with more direct approaches. 
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1 Introduction 



In this article we derive an equation for the cusp anomalous dimension for all angles and for 
all values of the 't Hooft coupling A = Oym-^ i'^ ^^^ planar limit of TV^ = 4 super Yang Mills. 
We obtain a system of non-linear integral equations of the form of a Thermodynamic Bethe 
Anstaz (TBA) system. The value of the cusp anomalous dimension can be obtained from 
a solution of the TBA system. This is also equal to the quark/anti-quark potential on the 
three sphere, see figure 1. 






(a) 

Figure 1: (a) A Wilson line with a cusp angle 0. (b) Under the plane to cylinder map the 
two half lines in (a) are mapped to a quark anti-quark pair sitting at two points on 5''^ at a 
relative angle of tt — 0. The quark anti-quark lines are extended along the time direction. 

The cusp anomalous dimension is associated with the logarithmic divergence arising from 
a Wilson loop with a cusped contour [1] 

(iy)~e"^^"=''^'^'^^'°^^, (1) 

where Ljr and euv are IR and UV cutoffs respectively. 

The locally supersymmetric Wilson loop in A/" = 4 super Yang Mills also includes a 
coupling to the scalar fields specified by a direction in the internal space n (with n^ = 1) 
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Instead of considering the same vector n on the two hnes that make the cusp, we can 
take two vectors n and n' . This introduces a second angle cos^ = n ■ n'. Thus we have 
the generahzed cusp anomalous dimension rcusp(0, 6*, A) [2]. Tcusp{(p,0) can be computed in 
terms of a solution of the TBA system of equations presented in this article. We can also 
consider the continuation = iip, where </) is a boost angle in Lorentzian signature. Before 
describing the computation, let us make some general remarks. 

1.1 Remarks on the cusp anomalous dimension 

Fcusp is related to a variety of physical observables: 

• It characterizes the IR divergences that arise when we scatter massive colored particles. 
Here ip is the boost angle between two external massive particle lines. For each consec- 
utive pair of lines in the color ordered diagram we get a factor of the form (1), where 
L/R is the IR cutoff and eyy is the given by the square of the sum of the momenta of the 
two consecutive particles. More explicitly, the angle is given by coshys = — ^7=|=?- This 

relation is general for any conformal gauge theory. See [3,4] and references therein. In 
TV = 4 super Yang Mills the massive particles can be obtained by setting some Higgs 
vevs to be non-zero $. Then the angle 6 is the angle between the Higgs vevs associated 
to consecutive massive particles [5]. 

• The IR divergences of massless particles are characterized by F^^p which is the co- 
efficient of the large ip behavior of the cusp anomalous dimension, Fcusp oc y^F^^p. 
F^gp was computed in the seminal paper [6]. Note that F^^p is also sometimes called 
the "cusp anomalous dimension" though it is a particular limit of the general, angle 
dependent "cusp anomalous dimension" defined in (1) . 

• By the plane to cylinder map this quantity is identical with the energy of a static quark 
and anti-quark sitting on a spatial three sphere at an angle n — (p. 

Fcusp(0,^)=V-(0,^). (3) 

See figure 1 . This potential depends on the angle as well as on the internal orienta- 
tions of the quark and anti-quark, which define the second angle 6. 

• In particular, in the small 6 = n — (p limit we get the same answer as the quark-anti- 



quark potential in flat space^ 

rcusp(0,A)~^^^^, when ,5 = vr-0^O, (4) 

where f (A) is the coefficient of the quark-anti-quark potential, V = ^^ ' , for a quark 
and an anti-quark at distance r in flat space and couplings to the Higgs fields which 
are rotated by a relative angle 6. 

• In the small (f) limit the cusp anomalous dimension goes as 0^ and one can define a 
Bremsstrahlung function B by 

rcusp~-(0'-^')5(A) 0,e«i. (5) 

This function B can be computed exactly using localization, see [8] and [9]. Here we 
will derive a set of integral equations that also determines B. In this way we can link 
the localization and integrability exact solutions. This function B is also related to a 
variety of observables, see [8, 9] for further discussion. 

Another motivation to study the cusp anomalous dimension is the study of amplitudes. 
Amplitudes are also functions of the angles between particles. Here we get a very simple 
function of one angle which has a structure very similar to amplitudes, since it is related to 
amplitudes of massive particles. Thus, obtaining exact results for this quantity is useful to 
learn about the general structure of the amplitude problem. 

1.2 Method 

The method to obtain the equation is a bit indirect and we need several preliminary results 
that are interesting in their own right. Just for orientation we will outline the main idea and 
method for its derivation. 

The method consists of the following steps 

• We ffist consider the problem of computing the spectrum of local operators on a Wilson 
line. We consider the particular case of operators with a large charge, i.e. operators 
containing a large number, L, of the complex scalar field Z insertions. These insertions 
create a BMN vacuum [10]. 



^This limit does not commute with the perturbative expansion in A. So (4) is correct if i5 ^ A. If we 
expand first in A and then take the 5 — >■ hmit we get a different answer due to IR divergences that arise in 
the naive perturbative expansion. These also arise in QCD, The origin of these logs are discussed in [5,7]. 



• In the large L limit the problem can be solved using an asymptotic Bethe Ansatz that 
involves the propagation of certain "magnons". These equations describe magnons 
moving on a long strip of length L with two boundaries associated to the Wilson loop 
on each of the "sides" of the operator, see [11]. The propagation of the magnons in 
the bulk is the usual one [12]. The new feature is the existence of a boundary. The 
magnons are reflected at the boundary and one needs the boundary reflection matrix. 
This is flxed in two steps. 

• We determine the matrix structure of the reflection matrix from group theory, as 
in [13-16]. This reflection matrix is such that it obeys the boundary Yang Baxter 
equation [17]. This is evidence that the boundary condition preserves integrabihty. 

• We derive a crossing equation for the reflection phase and we flnd a solution. 

• Doing a time/space flip, so that now we have the mirror theory between two boundary 
states separated by a mirror "time" L. See flgure 2. We can apply a symmetry generator 
that rotates one boundary relative to the other, so that we introduce the two angles. 

• We compute this overlap using TBA equations for any L, focusing on the ground state 
energy, which is extracted by taking the large T limit of the computation in flgure 
2. These boundary TBA equations can be derived following a method similar to the 
relativistic case [18]. 

• We set L = we get the cusp anomalous dimension. 

Let us discuss these steps in a bit more detail. First we should note that the exact inte- 
grabihty methods, as currently understood, work best to compute energies of states. Thus, 
we should phrase the computation of the cusp anomalous dimension as the computation of 
an energy. This is very simple. Under the usual plane to cylinder map, the cusp on the 
plane maps into two static quark and anti-quark lines on S^ x R. The quark and anti-quark 
lines are extended along the time direction, and they are separated by an angle vr — on 
the S"^, see flgure 1. The case = 0, which is the straight line in the plane, is mapped to a 
quark-anti-quark pair at opposite points on the sphere. If 6* = 0, this is a BPS conflguration 
and the cusp anomalous dimension vanishes exactly for all A. In fact, for 9 = ±0 we continue 
to have a BPS conflguration [19] and the cusp anomalous dimension continues to vanish. In 
general, the cusp anomalous dimension is the energy of this quark-anti-quark conflguration, 
as a function of the two angles, and 6. 

The conflguration with 6 = (p = preserves 16 supercharges which, together with the 
bosonic symmetries, give rise to a 05*^(4* |4) symmetry group. This is important to determine 
the boundary reflection matrix. 
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Let us begin by considering an apparently unrelated problem which is the problem of 
computing the anomalous dimension of operators inserted along a Wilson loop. First we 
consider a straight Wilson loop and we insert an operator at the point t = 0. For example, 
we can consider an insertion of a complex scalar field Z on the contour 

Pe'^--(^+**4)^(0)e'^o°°(^*+**4) ^ BiZ{0)Br . (6) 

These are operators that live on the loop and should not be confused with closed string 
operators. We denote these operators as BiZBj., where Bi^r stands for the usual path ordered 
exponentials of the gauge field. The operator considered above is BPS if Z is constructed 
out of scalars that do not appear in Bij. (6). To be definite, we consider Z = $5 + i^Q. 
We can similarly consider operators of the form BiZ^B^ which continue to be BPS. The 
straight Wilson loop is invariant under dilatations, so we can characterize the operators by 
their dimension under dilatations. These operators have dimension A = L. 

Determining the scaling dimension of operators of this type, but with more general inser- 
tions, is easier in the large L limit. Then, we can solve this problem by considering impurities 
propagating along a long chain of Z's. The impurities are the same as the ones that were 
used to solve the closed string problem in a similar regime [13,20]. The new aspect is that 
the impurities can be reflected from the boundaries at the end of the chain. This picture 
was discussed at the 1-loop order in the weak coupling limit in [11]. To proceed, we need to 
determine the boundary reflection matrix to all orders in the coupling. The matrix structure 
can be determined by group theory, as in [13]. The phase factor is more subtle. We write 
a crossing equation for it and we solve it following the strategy outlined in [21,22]. At this 
stage we have completely solved the problem for operators with large L. Up to corrections 
of order e~(^°'^*^*)''", we can find the energy of any open string state by solving the appropriate 
Asymptotic Bethe equations. 

After we have found the boundary reflection matrix we can then consider the possibility 
of rotating the half Wilson line that is associated with it. This rotation will simply act on 
the indices of the reflection matrix via a global transformation. Now we can consider states 
of the form BiZ^Br{6,4>), where we have rotated one of the sides of the Wilson line. This 
operator is no longer BPS but its energy is very small when L is very large, i.e. it has zero 
energy up to e~^^°^^^'>^ corrections. These are called Luscher (or wrapping) corrections. 

Before writing down the Thermodynamic Bethe Ansatz that describes the most general 
flnite L state, we will make some checks on the phase that has been obtained. As a non- 
trivial check one can get the flrst corrections to the ground state energy for large L. Namely, 
we are interested in the anomalous dimension of the operator of the form BiZ^Br{0, 0). This 
correction is given by a Luscher-type formula. This formula can be most simply understood 
by considering the problem in the mirror picture. Namely, we exchange space and time in 
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Figure 2: The BTBA trick. The same partition function can be viewed in two ways (7). In 
the open string channel it is a trace over all states in the open string Hilbert space. In this 
case Euclidean time runs along the T arrow. Alternatively we can view it as the propagation 
of a closed string along the L arrow. The closed string has length T and propagates over 
a Euclidean time L. The two boundary conditions, now lead to two boundary states that 
create the closed strings that propagate along the closed string channel. 

the open string picture. In other words, we have the equivalence 

^IX = Tropen[e-^^^^^^] = (^He'^^-- |i?.) , (7) 

where H'^'^^ is the open chain Hamiltonian on a strip of length L and -f/^cioscd is the closed 
chain Hamiltonian of the mirror theory on a circle of size T. So now we have a closed 
string exchanged between two boundary states. The analytic continuation of the boundary 
reflection matrix gives us the probability of emitting a pair of particles from the boundary 
state. It turns out that this continued reflection matrix has a pole at zero mirror momentum 
which implies that we can create single particles [17]. The coefficient of the pole in the 
reflection matrix at zero mirror momentum determines the prefactor of the Luscher correction 
[23]. We compute this at strong coupling and we find agreement with a direct string theory 
computation. Furthermore, the leading order correction at weak coupling, going like g"^, also 
comes from this Luscher type term. In this way we match the leading corrections at weak 
and strong coupling. This constitutes a test of the boundary refiection matrix. In particular 
the very existence of the pole at zero mirror momentum is due to the phase factor of the 
matrix, which we derived by solving the crossing equation. 

Finally, one can write down a Thermodynamic Bethe Ansatz equation that describes the 
finite L situation. This follows the standard route for getting the energies of states of an 
integrable field theory with a boundary. The derivation of these equations is very similar to 
the derivation of the equations for closed string states. The new element is that instead of a 
thermodynamic partition function we have the overlap between two boundary states, as in 
(7). The derivation of TBA equations for integrable systems with a boundary was considered 
in [18]. The boundary states are given in terms of the analytic continuation of the boundary 
reflection matrix. The TBA system of equations arises from evaluating this exact overlap 



between the two boundary states in an approximate way by giving the densities. Most of 
the TBA equations come from the entropy terms, which are the same in our case. Thus the 
boundary TBA equations are very similar in structure to the closed ones. We obtain 

log Ya = log(4fi::4) - 2LEm,A + Kab * log(l + Yb) . (8) 

The cusp anomalous dimension, or quark/anti-quark potential is given schematically by 

oo 
"" ^ 

Here Em,A and qa are the energies and momenta of the excitations in the mirror theory. 
The equations will be given below in their full detail, (62)- (66). The information about the 
boundary is contained in ka which comes from the reflection phase of the theory and depends 
on the boundary state. 

Note: We were informed that similar ideas were pursued in [24]. 

2 Spectrum of operators on a Wilson line 

Let us first discuss the symmetries preserved by a straight Wilson line. Let us start with the 
bosonic symmetries. It preserves an SL{2) x SU{2) x 50(5) symmetry group. The SO{b) is 
the subset of 5*0(6) that leaves $^ invariant, where ^^ is the scalar that couples to the Wilson 
line. The SU{2) factor corresponds to the spatial rotations around the loop. The SL{2) 
factor contains time translations, dilatations and special conformal transformations along 
the time direction. In addition, we preserve half of the supercharges. The full supergroup is 
OSp{A*\A). The star means it is the real form of SO{A) such that SO{A*) ~ SL{2) x SU{2). 

Now we can consider the insertion of an operator of the form Z^ on the Wilson loop, we 
can denote this as BiZ^Bf Here we choose Z to be Z = $^ + i$^. 

The operator Z inserted at the origin preserves an SU {2\2)'^ subgroup of the full symmetry 
group of the theory. The Wilson loop, together with the Z insertions at the origin preserve an 
SU{2\2)j:, subgroup of all the symmetry groups we mentioned. This is a diagonal combination 
of the two SU{2\2) factors preserved by Z. This common preserved symmetry is very useful 
for analyzing this problem. These operators are BPS, and they have protected anomalous 
dimension, £ = A — J^q = 0. 

Note that on S^ we have a flux tube that goes between the quark and the anti-quark. 
These operators inserted on the Wilson loop are mapped to to various excitations of the flux 
tube. 



2.1 The boundary reflection matrix 

Recall that the bulk excitations are in a fundamental representation of each of the two iM(2|2) 
factors of the sm(2|2)^ symmetry of the Z-vacuum. The tilde means that we are considering 
the momentum dependent central extensions discussed in [13,25]. In other words, we can 
think of them as particles with two indices \E'^ ^, where A labels the fundamental of the first 
5u(2|2) and B labels the fundamental of the second su{2\2) factor of the sm(2|2)^ symmetry of 
the infinite chain. This central extension determines the dispersion relation for the excitations 
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The scattering matrix between two particles has the form S'^^J^ = SqS'^^S'^^ [13]. Namely, 
it is the product of a phase factor Sq and two identical matrices, one for each sm(2|2) 
factor. These matrices are fixed (up to an overall factor) by the ?u(2|2) symmetry of the 
theory [13,25]. These matrices depend on the two momenta, pi and p2, of the scattered 
variables. The phase factor 50(^1,^2) was guessed in [6,26], and a nice derivation was given 
in [21,22]. 
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Figure 3: Unfolding of R{p) into S{p, —p). There is a non-trivial map between dotted and 
checked indices. See appendix A for details. 

In our problem we need to fix a reflection matrix of the form R^^{p). Let us consider 
first the refiection from the right boundary, see figure 3. This matrix depends on only one 
momentum p, the momentum of the incident magnon. The boundary is invariant under an 
sm(2|2)d symmetry group, which is diagonally embedded in the 5u(2|2)2 symmetry group of 
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the bulk of the spin chain (see appendix A). A similar problem was studied in [16] and the 
matrix part of the reflection is the same. Thus the symmetries constraining the reflection 
matrix are exactly the same as those constraining the bulk scattering matrix for each of 
the sti(2|2) factors. From this argument we expect that the matrix structure should be 
completely fixed. In fact, the matrix structure should be essentially the same as what we 
encounter in the matrix S^j^{p, —p), or i?^^(p) oc S^'^{p, —p). One is tempted to say that 
the scattering phase factor would be 5*0(^5, —p)- However, this is not fixed by the symmetries, 
and will not be true as we discuss below. In the presence of a boundary, we can do a kind 
of "unfolding" of the spin chain. Here each bulk magnon is viewed as a pair of magnons of 
5^(2 1 2)d, one with momentum p to the left of the boundary and one with momentum —p to 
the right of the boundary. See figure 3 . 

This completely solves the problem of fixing the matrix structure of the reflection matrix. 
The full reflection matrix, in complete detail, is given in appendix A. One can also check 
that it obeys the boundary Yang Baxter equation. But this is clear from the "unfolded" 
picture in terms of a single chain. We should emphasize that we have assumed that there 
are no boundary degrees of freedom. We do not see any evidence of any boundary degrees 
of freedom at either weak or strong coupling, so this is a reasonable assumption. 

Before we determine the phase, let us make a side remark. There is a variety of problems 
that give rise to a spin chain with boundaries and preserve the same symmetries, 05^(4*14). 
We can consider an open string ending on a D5 brane that wraps AdS^ x S"^, or AdS2 x S^. 
In fact, there is a whole family of BPS branes of this kind that arises by adding flux for the 
U{1) gauge field on the brane worldvolume on the S"^ or Ad 82- In fact, in the limit of large 
electric flux on the AdS2 x S"^ brane we get a boundary condition like the Wilson loop one. 
In fact the AdS2 x S'^ branes can be interpreted as Wilson loops in the k-iold antisymmetric 
representation of U{N) [27]. In all these cases one can choose the BMN vacuum (or choose 
the field Z) in such a way that we preserve the 'su{2\2)d of the spin chain. Therefore, we 
would get the same matrix structure for the reflection matrix, again assuming that there are 
no boundary degrees of freedom. However, they would differ in the choice of a phase factor. 
Below we get a phase factor which has all the right properties to correspond to the one of 
the Wilson loop. It would be interesting to fix the phase factor also for these other cases, 
but we leave this to the future. 

In order to fix the phase factor we write a crossing equation. We derive this by writing 
the identity state of [13], scattering it through the boundary and demanding that the full 
phase is equal to one. Denoting the phase factor as Rq, defined more precisely in appendix 
A, we obtain the crossing equation 

Ro{p)Ro{p) = a{p,-p)\ (13) 
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where the bar indicates the action of the crossing transformation. Here (y{pi-,P2) is the 
bulk dressing phase, discussed in [6,22]. We are going to p along the the same contour in 
momentum space that we choose in the formulation of the bulk crossing equation. 

In addition, we also should impose the unitarity condition 

R,{p)R,{-p) = l. (14) 

We now write the ansatz 

1 A 



RJp) = -—-^, I ^ 1 . (15) 

Here a is the bulk dressing phase. This would be our naive choice for a phase factor. The 
explicit factors of x^ have been chosen only to simplify the final formula. We have an 
unknown factor crsip)- Now (13) becomes 

x~ + :^ 



aB{p)aB{p) = ::^y^ . (16) 



We can now solve this equation using the method proposed in [21,22]. We give the details 
in appendix B. We obtain 

aB = e*^(^^)-*^(^"^ (17) 

• . N .^z N I dz 1 ^ rsinh[27r^(z + ^)]^ , , ^ 

tx(x) = t^(x)= (b log<^ -^^-^T^>, x>l. 18 

\z\=l 

This expression is valid when |x| > 1. The value for x in other regions is given by analytic 
continuation. We have also introduced the function $(x) which is given by the integral for 
all values of x. When |a;| < 1 these two functions differ by 

rsinh[27ro(x+ -)] 1 
tx{x) = z^x) + \og\ ^ ;\ ^■''^ , x<l 19 

The ambiguities in the choice of branch cuts for the logarithm cancel out when we compute 
aB in (17). Note that xi^) = x{~x)- 

So far, we have found a particular solution of the boundary crossing equation. Still, the 
true dressing phase might require the inclusion of further CDD factors. In order to make a 
conjecture for the exact boundary dressing phase, we need to compare against some explicit 
computations. 

Before doing so, let us observe that, given (Tb{p), we can define an infinite family of 
solutions by taking 



(P) = ( ^T^ ) i^^(P)] ^ • (20) 
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By computing the dressing phase in the physical regime we will be able to show that s = 
is the solution we want. 

The proposal for the phase factor for the reflection matrix, given in (15), (17) is one of 
the important results of this paper. We will perform various checks on its validity. 

2.2 Checks of the boundary reflection phase in the physical region 



P = TV 
P / \ \ -P 
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Figure 4: Computation of the reflection phase at strong coupling. We have a soliton at the 
boundary, which is at rest at a = 0. There is also an image soliton coming from the right. 
Then the soliton with momentum p scatters through the soliton at rest and the one with 
momentum —p, leading to a certain time delay. From the time delay we can compute the 
derivative of the reflection phase with respect to the energy. 

Let us describe how to compute the boundary dressing phase at strong coupling. We 
have to consider the open string solution that corresponds to a 1/2 BPS Wilson line carrying 
a large J^q charge given in [11]. This solution describes the transition from the boundary 
Wilson line to an inflnite BMN vacuum. It is convenient to understand this solution in 
the conformal gauge, when we set the stress tensor on the S^ equal to one, and the stress 
tensor of the AdS^ to minus one. In these variables, the problem only involves an AdS2 x S'^ 
subspace and we can perform the Pohlmeyer reduction in each factor. The S"^ part gives 
rise to a sine gordon theory and the solution is just half of a soliton at rest. More precisely, 
the center of the sine gordon soliton sits at the boundary. In the AdS2 part we have a sinh 
gordon theory, and the solution is a sinh-gordon "soliton". This is a singular solution which 
is the direct analytic continuation of the sine gordon soliton. The singularity reflects the fact 
that the string goes to the AdS boundary. If we compute the energy, there is a divergent 
part and the flnite part is zero. The setup is explained in more detail in [11]. The fact that 
the flnite part of the energy is zero is consistent with the absence of a boundary impurity 
transforming non-trivially under sm(2|2). A bulk magnon is a sine gordon soliton, and leaves 
the AdS part of the solution unperturbed. In the presence of a boundary, we need to put 
also the "image" of this soliton and the conflguration looks as in flgure 4. The reflection 
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involves the scattering of the sohton with the image sohton as well as the scattering with 
the soliton at rest. These soliton scattering phases were computed in [14,28]. So the strong 
coupling limit of the right boundary scattering phase -Ro(p) = e*'^^'-^-' is given by 

/ 1 _ sin 2 \ 
Sr{p) = -8g cos 2 log cos f - Ag cos f log ( ^^—^ j ■ (21) 

The first term in (21) is exactly what one gets from the strong coupling limit of the factor 
l/cx^p, —p) [28]. We will see that the second term corresponds to asip)^^- At strong coupling 
can expand (18) as 



ix{x) ~ ig 
which leads, for physical excitations, to 



1, 1 , (x + i) 

x+-)-log) ^ 

X 2t [x — t) 



+ 0(1) , (22) 



log an {p) = Ag cos | log ( —^ ] ■ (23) 

1 + sm - 



1 , , . . . , / 1 - sin 2 ' 

2 . 

This indicates that we must pick the case s = from the family of solutions (20). 

Finally, let us discuss the behavior at weak coupling. The bulk dressing phase a has its 
first contribution at order g^ , leading to four loop corrections to anomalous dimensions. On 
the other hand, the boundary dressing phase, as, receives its first contribution at g^, so that 
it will start modifying anomalous dimensions of operators inserted on the Wilson loop (dual 
to open string states) at three loops. 

2.3 Reflection matrix for a Wilson line at general angles 

We will need the boundary reflection matrix for a Wilson line sitting at general angles, and 
9. In particular, we want the left and right boundaries of the open chain to be rotated by 
relative angles. We can obtain the boundary state of the Wilson line at a different position 
on the S*^, or the S"^, by applying a symmetry transformation on B,.. This should be a 
symmetry that is broken by Br- So for example, we can apply an SU{2)l rotation on the 
S^ which is in one of the SU{2) factors in the S0{4) group of rotations of the 3-sphere. If 
we apply an SU{2)l rotation with an angle 20, we will get that the quark is rotated by an 
angle on the S^, away from the south pole. See figure 1(b). Note that the SU{2)l we are 
considering is a symmetry of the Z vacuum. We have a similar feature on the S^ . We can 
also apply a rotation in an SU{2)li factor inside S'0(4) C 5*0(6) (this 5*0(4) leaves the Z 
vacuum invariant). The reflection matrix is very easy to obtain. We pick these two SU{2) 
generators so that they sit in the bosonic part of one of the sm(2|2) factors of the sm(2|2)^ 
symmetry of the bulk. Then they will simply introduce some phases of the form e*"^ or e*^ 
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when a state is reflected from tlie boundary and its SU{2) quantum number clianges. Tlie 
SU{2)l or SU{2)li quantum numbers can cliange because tliey are not symmetries of the 
boundary state. More explicitly, the reflection matrix from a boundary state at angles 0, 9 
is given by 

i?Jj(^, 0) = (m-i)£m5i?g|(0, 0) , with m = diag(e*^ e''^, e'^, e"*^) . (24) 

Note that the matrix m acts only on the undotted indices since we did a rotation inside only 
one of the Hi(2|2) factors. 

2.4 Luscher computations and checks in the mirror region 

In this subsection we start considering the problem with two boundaries. In other words the 
operator BiZ^Br{(f),9). Here (p, 6 are the relative orientations of the two boundaries. On 
the plane, this corresponds to a cusp, plus an operator of the form Z^ at the tip. In the 
limit L ^ 1 we get the naive superposition of the two boundaries and the energy of the state 
is zero {S = A — L = 0), regardless of the orientation of the two boundaries. The leading 
correction is of the form e-(constant)L 'jj^ggg corrections come from the exchange of particles 
along the "mirror" channel. The boundary sources particles, which then travel to the other 
boundary. These corrections sometimes go under the name of Luscher corrections. Of course 
the familiar Yukawa potential is a simple example where the leading correction comes from 
the exchange of a single massive particle. 

In order to derive the precise correction formula it is convenient to describe in more 
detail the mirror theory. In the bulk of the worldsheet the mirror theory was discussed in 
various papers, see [29] for example. This theory is defined by exchanging the space and time 
directions of the spin chain we have been considering so far. Thus, instead of (11) we define 
q = ie and Em = ip, and use the same formulas as in (11). Here q is the mirror momentum 
and Em is the mirror energy. In order for these to be real we will need to pick a solution of 
(10) with |x'''| > 1 and |x~| < 1. From the expression for q, we can write 

,[±^] = _^i /i + ^r;L^±i)(g + ,a), (25) 




Em = 2arcsinh-^^ . (26) 

Here z^ just denote the values of x^ in the mirror region. We have also written the dispersion 
relation in the mirror region, for an arbitrary bound state. The elementary mirror magnon 
has a = 1. 

When we have a boundary, this time/space flip turns the boundary into a boundary 
state, see figure 2. Then a suitable analytic continuation of the boundary reflection matrix 
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characterizes the boundary state. The boundary state creates a supersposition of many 
particles. The total mirror momentum should be zero since it is translational invariant. So, 
schematically the state has the form 

oo 

IB) = 10) + / I^I<''-''%)-l,,A<ss\0} + ■ ■ ■ (27) 



with 

irMBB(^) = [i?-i(z+,^-)];^jc:^^'^^, (28) 

where we put the mirror values (25). Here R is the right reflection matrix, with z^ contin- 
ued to the mirror region (25). This amounts to an analytic continuation of the reflection 
matrix. Here cl ■ is the creation operator of a magnon with momentum q. C is a charge 
conjugation matrix. In the case of a relativistic model with a single particle (28) reads 
K{6) = 1/R{6 — z|), [17]. The formula (28) can be obtained by performing a 7r/2 rotation 
of the boundary condition. Due to the independence of reflection events from a boundary, 
we can exponentiate (28) to get the full boundary state [17, 18]. Similarly, we can form a 
future boundary state. This is a boundary state that annihilates the particles. It is given by 

oo 

{B\ = (0| + (0| I ^K^A,BBiqK^<^-f + ■■■ (29) 



with 







R- 



1 1^ ^^ 



z z^ 
1 



^ddaA ■ (30) 

BB 



In the relativist case (30) would be K{6) = j^r^T^^^y 

When L is very large the leading L-dependent contribution comes from the exchange of 
this pair of particles and we can write the corresponding contribution as 

oo 

^^ = '1 ^e-2^^-('')t(g) , t{q) = TT[K{q)K{q)] . (31) 



This formula is correct whenever the integral is finite. 

In our case, the phase factor cxb has a pole at g = 0. In the physical region (Tb{p) was 
perfectly finite. This pole in the mirror region is crucial for obtaining the correct answer. But 
first we need to generalize (31) to the situation when we have a pole at g = 0. The physical 
interpretation of this pole at g = is that the boundary state is sourcing single particles states 
in the mirror theory [17]. For a similar case in the AdS/CFT context see [30]. Obviously 
such source has to contain only zero momentum particles. 
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A careful analysis leads to the formula [23] 

oo 

^ ~ - / ^ log {1 + e-^^^™('^)Tr[ir(g)i?(g)]} ~ -^e-'^''-^''^ ^ q^i:v[K{q)K{q)]\,=o ■ (32) 



In the last equality we extracted the leading term in the integral, which comes only from the 
coefficient of the pole. Notice that the L dependence is precisely what we expect from the 
exchange of a single particle. We should sum over all the particles that can be exchanged. 
The mirror theory contains bound states indexed by an integer a = 1, 2, ■ ■ ■ , and we should 
sum over them. 

In appendix D we show that we can evaluate t{q) for a fundamental mirror particle and 
we obtain 

t{q) = as{z\z-)aB{-^,-^^ (j:) ^ (Tr[(-lf ])% (33) 

where the trace is over the four states of a single 'su(2\2) magnon. Let us now give a simple 
explanation for this formula, for more details see appendix D. We can write the reflection 
matrices that appear in K and K (28) (30) in terms of bulk S-matrices for the unfolded 
theory, namely in terms of bulk 5* matrices for a single ?u(2|2) factor. The matrix in K is 
essentially S{—p,p) and the one in K is S{p, —p). When we multiply these matrices we can 
use the bulk crossing equation to get the identity. Here we should use the full bulk matrix, 
including the bulk a factor. This is the reason that the bulk a factor disappears from the 
final formula (33), but the boundary one remains. The factor of z~ / z^ arises from the factor 
in parenthesis in (15). Finally, the (—1)'^ is related to the fact that we have fermions. Here 
F is the fermion number. When we perform the TEA trick, we get periodic fermions in 
Euclidean time if we started with periodic fermions in the spatial directions. Of course a 
periodic fermion in Euclidean time is the same as the trace with a (—1)^ inserted. The 
operations that lead to (33) can be understood graphically as in figure 5. 

Of course, for a fundamental magnon Tr[(— 1)^] = 0. This is good, since it is saying that 
the correction vanishes in the BPS situation. If we rotate one boundary relative to the other 
then we need to perform the replacement 

Tr[(-1)^] — > Tr[(-l)^m] = -2(cos0 - cos^) . (34) 

where m is given in (24). Again, we see that it vanishes in the BPS case = ±0. 

To write down the full Luscher formula we need to compute t{q) also for the bound states 
of the mirror theory. One can first use the standard fusion procedure to get the bound state 
refiection matrix. Then one can use the same argument as above to eliminate the bulk S 
matrices, as in figure 5. The final formula is 
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(a) 



(b) 



Figure 5: (a) We have a strip with pairs of particles being exchanged. The two colors 
represent the two types of indices. In (b) we unfolded this into a cylinder computation. The 
K matrices became 5* matrices for a single su(2|2). (c) Using crossing we have moved the 
lines. The red circles indicates the action of the matrix m. 



ta{q) 



aB{z^^''\z^-^^)aB 



y[+a] 



(Tr[(-l)^mj)^ 



Tr[(-l)^mJ = (-l)"2(cos0-cos( 



smac 



smi 



(35) 
(36) 



where now the trace is over all the states of a magnon boundstate in a single copy of iM(2|2), 
see equation (167) in appendix D. As anticipated, an important property of cr^ is that it has 
a pole at g = 0. More precisely the combination of (Jb in (35) becomes 

1 N 



3ix(z[+''l)-Jx(2'"''')+«x(l/^'~''')-Jx(l/2'+''') 



2'Kg{z^-''\ 



2'Kg{z^+^^ + ^\ 



sinh[27r(?(z[-«] + -i^)] sinh[27r(?(zM + ^)\ 



xe 



i(<I>(2[+'^l)-*(2[-°l) + #(l/zl-"l)-${l/z' + ''l)) 



(37) 



Here we used that z is in the mirror kinematics and we used (19) to evaluate x{^) when 
|a;| < 1. We have also used that x{~^) = x(^)- Each of the sinh factors leads to a pole at 
g = 0. Namely, using (25) we get 



2vr^(^'^"l + ^) = i^vra + Tiq\jl + ^ + 0{q') 
for small q. We then can write the pole part of (37) as 

1 4 



q^ (a^ + 16^^^ 



-F{a,gY + 0{l), 



with 



F{a,gf 



^i(<I>{2[+'"l)-$(2[-''l) + <J>(l/2[-''l)-*(l/2[+''l)) 



g=0 



(3^ 



(39) 
, (40) 



where the last factor is evaluated at g = 0. 

Then we find the coefficient of the double pole of t as 

{cos (j) — cos Oy . 2/ ,N o,^ / — a + a/oM-IO^ 



2 



\imqHa{q)=4'- ^. ^ sin^(a0) ,_^.^,. [ j====^ ] F{a,g)\ (41) 

9^0 sm^ (f) [a^ + I6g^) \ a+ J a^ + \%g^ 



The factor in parenthesis is [z^ "l/^t"')^. Finally, inserting this into the expression for the 
energy (32), we find 

.^ (cos0 — cos6') ^, ^ , _ / — 1 + \/l + \&g^ ja^ \ . / ,x a i^, \ 

A^~-^ ^-^ — ; ^> -1 , -\ sm a0^^= F{a,g) ■ 

sin0 4^' \^ 1 + v^l + 16^Va2 ^ ^ v^TTT^V^ 

(42) 
The factor in parenthesis is just ^-^rn.{a){L^-\) ^ representing the exchange of a bound state in 
the mirror channel. The sign (—1)" is a bit subtle and has to do with the correct sign we 
should pick for the square root in (32). The correct sign is easier to understand for an angle 
of the form = vr — 5, for small 5. In this case we have a quark antiquark configuration and 
it is clear that we should get a negative contribution to the energy. In fact, we can think of 
the overlap of the two boundary states as computing a kind of norm or inner product. We 
see that in terms of b the expression has the expected sign. In other words, for small b we 
get the positive sign of the square root in (32). Of course, once we get the expression for 
small (5 we can write it in terms of 0, or even analytically continue = —i^p. 

2.4.1 Leading Luscher correction at weak coupling 

The expression (42) gives the leading Luscher correction at all values of the coupling for 
large L. Let us now examine it at weak coupling. Then the factor in parenthesis in (42) is 
of order g^ . So, at leading order, we get a term of the form g^^'^^ . This has the interesting 
implication that this leading "wrapping" correction appears at L + 1 loops. In particular for 
L = 0, the one loop contribution comes from such a term!. In fact, expanding (42) to leading 
order in g^ ^ and setting L = 0, we can set F = 1 to this order and obtain 



sm.(i) ^-^ a 

a=l 

r, (cos d) — COS 6*) , , , 

= 2^2 1 V 1^^ 44 

sm0 

which coincides exactly with the leading 1-loop contribution to rcusp(0, ^) computed in [2]. 
We can also do the computation of the leading order term for any L, we get 

E = _,2.2.(cos0-cosg)(-l)^(4vr)|-^^^^^^^^ /vr-0\ ^ ^^^,,,,^ ^^^^ 



sincj) (1 + 2L)! "^^"^ \ 2n 

19 



where -Bn(x) is the Bernoulh polynomial, which is a polynomial of degree 2L + 1. In [5] a 
particular class of diagrams was identified which produced the same expression. 

2.4.2 Leading Luscher correction at strong coupling 

We can also compute the leading large L correction at strong coupling. We simply evaluate 
the large g limit of (42). First we note that 

L+l 

^ e-*^ = e-^^™('^=°) . (46) 



This implies that to leading order in e~^ we only need to consider the case a = 1. The 
expansion of the function F is done in appendix C.l eqn. (160). Putting everything together 
we find that the leading strong coupling correction goes as 

£ = (cos - cos 6*) ^e~ ^ . (47) 

This agrees precisely with the result computed directly from classical string theory in ap- 
pendix C, see (152). This constitutes a nontrivial check of the reflection phase. Notice, the 
funny factor of e~^ which is correctly matched. 



3 The open Asymptotic Bethe Ansatz equations 

We will now write down the asymptotic Bethe ansatz (ABA) equations that describe the 
spectrum of operators with large L inserted on the Wilson loop. These give rise to a spin 
chain with two boundaries, which are separated by a large distance L. Moreover, the ABA 
equations are used to derive the BTBA system by embedding them into the closed equations, 
as we do in appendix E.l. 

In order to obtain the ABA equations we have to diagonalize the way the bulk and bound- 
ary scattering matrices act. This can be done by formulating a nested Bethe ansatz, which 
defines impurities at different levels of nesting. Here we just sketch the computation, which 
is a straightforward generalization of the case with periodic boundary conditions studied 
in [13]. 

Consider an asymptotic state with A^^ bulk magnons, or level I excitations, on the half- 
line with a right boundary. We will introduce a second boundary and relative angles later, 
when writing down the Bethe equations. In particular, we can consider a state whose level I 
impurities all carry the same SU{2\2)d index'^. Say, for example, in the unfolded notation, 

|*3(pi) ■ ■ ■ ^3{pm)^3{-Pm) ■ ■ ■ ^3(-Pi)) = |0)" , (48) 



^The choice of index is arbitrary. 
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which is regarded as the level II vacuum state. Of course, we could also consider states where 
A^^^ out of the iV^ level I impurities have different indices. Those should be understood as N^^ 
impurities in the level II vacuum state. In total, such states will contain A^^ level I impurities 
and A^^^ level II impurities. In general, we have I'^aiivi) ■ ■ ■ ^^aj^uiym^))^^ for o-fc = l, 2, where 
i/k are auxiliary parameters associated with the level II impurities. 

Similarly, a third level of nesting can be defined. If all the level II excitations carry the 
same index, for instance |\l/i(yi) ■ ■ ■ ^i(y7viO)^^ = |0)^^^ "^^ can define a level III vacuum 
state. Then, magnons \l/2 will be treated as level III impurities propagating in |0)^^^. For the 
kind of SU {2\2) spin chain we are considering, this level III is the final level of nesting^. 

Then, to formulate a coordinate Bethe ansatz, bulk and boundary scattering factors 
among excitations of different levels have to be introduced to write the nested wavefunc- 
tions. Those can be determined by imposing certain compatibility conditions. Namely, that 
the action of the bulk and boundary scattering matrices on wavefunctions with higher level 
impurities just pulls out the same factor as when acting on the level II vacuum state. Nat- 
urally, the bulk scattering factors are exactly the same as the ones obtained in the periodic 
case [13], 

S''\xf,xt) = -Soipi,P2), (49) 

5^'"(x^y) = 1/S'''\y,x^) = -^^, (50) 

S'''''\w, y) = ^ ' I' '! = ^ , (51) 

W — V — TT 

S'''^''\w^,w,) = —I, (52) 

Wi -W2 + 



w - 


-y- 
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-y- 
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Wi- 


-W2 
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9 





9 



where 






is the bulk dressing factor and v = y + -. All other bulk scattering factors are trivial. 

The reflection factors can be derived in the same way. The level II vacuum (48) containing 
A^^ magnons becomes a lattice with 2N'^ sites. Consider a single level II impurity propagating 
in this vacuum from left, i.e. propagating along the left (undotted) indices of the bulk 
magnons. Undotted and dotted indices can only mix by the reflection of the rightmost bulk 
magnon. That could make us think there exists a defect in the middle of the level II vacuum 
lattice which separates the 3 and 3 indices of the rightmost level I impurity. In principle, 
the level II impurity could be reflected and transmitted across such defect, see figure 6. 

^\I>4 arc not considered as elementary but as double excitations. 
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Figure 6: Propagation of a single level II impurity across the defect. 



However, and because the boundary scattering matrix R{p) oc S{p, —p), the compatibility 
condition we obtain from the reflection of the rightmost level I impurity is analogous to the 
ones we obtain from the scattering of two level I impurities. In this way, the compatibility 
conditions imply that level II impurities are purely transmitted. In other words R^^ = and 
R^^ = 1. Analogously, the reflection of level III impurities is determined. In summary, we 
have 

R\x^) = Ro{p), R'\y) = l, R''\w) = l, (54) 

where Ro{p) is the boundary phase factor (15). 

Let us now put the system in a finite strip by introducing another boundary. We will 
then have certain quantization conditions on the rapidities for all kind of excitations, namely 
the Bethe ansatz equations. We will introduce the left boundary with relative angles with 
respect to the right one, by using the rotation discussed in section 2.3. To understand how 
this rotation affects the factors R^, R^^ and R^^^ it is enough to consider the action of m, 
defined in (24), on the following key components of the reflection matrix 



p33 , p33 


=^ 


R^^R\ 


Ri ^ e^'-^R^ , 


=^ 


^11 ^ g.e-.0^ii 


21 ' ^ 21 ' 


=^ 


Rin ^ e-2.e^iii 



(55) 



Let us finally write down the nested Bethe ansatz equations. They are obtained by 
picking an impurity of any level of nesting and moving it through all the other impurities 
twice and reflecting it from both boundaries as it is shown in the left picture of figure 7. 

If we go to the unfolded picture what we have is periodic chain of length 2L, where for 
every level I excitation of momentum pk there exists a mirrored one of momentum —pk, figure 
7. Such duplication does not occur for higher levels of nesting, for which the excitations do 
not necessarily come in pairs. When moving around the level I excitations to derive the 
Bethe equations, we have to recall that their duplication is an artifact of the unfolding. 
Every pair represents a single magnon in the original picture. When we move the original 
magnon, it looks like moving the pair simultaneously in the unfolded picture. Then, for level 
I impurities we pick up the factors that correspond to simultaneously moving around the 
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Pi P2 P3 P Pi -Pi -P -P3 -P2 -Pl 

(b) 

Figure 7: Bethe equation for the open chain, (a) The original picture with boundaries. 
The particle goes to one boundary, then the other, and finally back to the origina position. 
(b) The unfolded picture. We have a closed circle. The leftmost solid line is identified the 
rightmost one. The motion that leads to the Bethe equations involves moving the magnon 
with momentum p around the closed circle and at the same time we also move its partner 
which has momentum —p around the circle in the opposite direction. 



pair with momentum pk and —pk in opposite directions. For level II impurities, we have to 
collect the factors corresponding to going through all the level I pairs and all the level III 
impurities (scattering between level II particles is trivial). Finally, for level III impurities we 
get the factors of going through all the level II impurities and all the other level III impurities 
(scattering between level III and level I particles is trivial). The resulting set of open Bethe 
ansatz equations is the following 
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^b{p 



kr<y{pk, 
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N' 
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(56) 



(57) 
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Eq. (56) can be re-written, including I = k in the second product, as 

2L _ , 1 ^ Af" 



(+X 2_L _ , 1 ^ JV" 



£1 (4 - ^r)(l - :;^) (^^ + 4)(1 + :;^) 



As usual, the energy is given by 

£ = Y,<Vk) (60) 

fc=i 

4 The boundary TBA equations 

The Bethe equations (56)-(58) presented in the previous section are the correct description 
of the spectrum for large chains, L ^ 1. As L becomes small, wrapping effects come into 
play and the Bethe equations are no longer valid. Moreover, in this paper, we are mainly 
interested in L = 0. A description of the spectrum that is valid for any L is the Boundary 
Thermodynamic Bethe Ansatz (BTBA) equations. These are a set of integral equations 
that govern the dynamics in the mirror channel. That is, the dynamics of excitations after 
exchanging the two dimensional space and time directions [17,31], see figure 2. The TBA 
equations can be derived from the knowledge of the spectrum of states and bound states in 
the mirror channel. This spectrum was derived in [32]. The derivation of the TBA equations 
then follows the standard route given in [31,33-36]. In the case that we have a boundary we 
can follow essentially the same route. We use the boundary state defined in section 2, and 
the untangling of boundary reflection matrices described in figure 5. Then we get a TBA 
which looks very similar to what we would obtain for a closed chain of twice the length L, 
except for the fact that for each particle of momentum q we get one of momentum — g, since 
the boundary state creates such a pair of particles. The consequence of this is that the Y 
functions obey a reflection property 

Ya,s{u) = YaM-u) (61) 

The set of Ya^s functions is the same as the one we have for the closed string problem 
[34-36]. However, due to (61) we can restrict our attention to the ones with s > 0. The 
boundary data appears as chemical potentials which depend on the angles, 6, 0, as well as a 
u dependent chemical potential given by the boundary dressing phase as- The precise form 
of the equations is derived in appendix (E). 
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Figure 8: (a) Set of Fa,s functions for the closed string problem. Here we have the same set 
but the additional condition (61) implies that we can restrict to the set in (b). 



Let us summarize the final equations 

^1,1 r^ 1 1 + Yim 1 + Ym 1 
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^2,2 
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Km-l * log 
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2,2 



where we used the conventions of [37,38] for the kernels and integration contours'"*. We have 



also defined the barred F's as Y 



(here) 



IK 



(there) 



(see appendix E.3 for a summary). Here, 



the momentum carrying Y^ o functions are defined as symmetric functions Yafl{—u) = Yafi{u) 
and * /(f) = [*f{v) + */(— f)]/2 is a symmetric convolution^. There are implicit sums over 



^ The convolutions of terms depending on Yi i or 1^2.2 are over a finite range |m| < 2g. We use * as a 
reminder of that. 

^For the ground state, we expect all functions to be symmetric, Ya.s{u) — Ya.s{—u). But for excited states 
(61) only requires the Ya,o functions to be symmetric. The equation for excited states could in principle be 



obtained by analytic continuation from these equations [37-40]. 
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one of the indices of the kernels^. The bold face Y's represent the asymptotic large L 
solution. This is the solution we obtain when the convolutions with the momentum carrying 
YaflS are dropped. These asymptotic solutions are the only place where the angles and the 
boundary dressing phase enter. They are given by 

^ _cos^ - ^ sin^^ 

''' COS0 ' '" sin[(s + 1)0] sin[(s - 1)^] ^ ' 

Y22 = -^^, Y,i = ""^^ (68) 

COS0 ' °' sin[(a + 1)0] sin[(a — 1)0] 



Y„n = 4 



— 71 — 7 — , , r I u , , 1 (cos d) — COS y) ^ — 



where x is the function defining the boundary dressing phase (17). Notice that the length L 
appears only in (69). Here z^^°^ are the solutions of 



^ = M^'''"' + -r^)-^^ = M^'~"' + -r^)+^^ = :^Wi + -ff^ (70) 



in the mirror region with \z^'^'^\ > 1 and Iz^'""'] < 1. 

Once we solve this system of equations, we can compute the ground state energy as 



^ = -E/|?Ml + X.,o), (71) 



a=l Q 



where q is the mirror momentum of each magnon bound state 

1 



q = 9 






(72) 



4.1 Recovering the Luscher result 

As a simple check of these equations let us rederive the results of section 2.4. In the large 

/ r_ , \ 2L + 2 

L limit we see that the factor ( ^pp^ j = Q-^rn2{L+i) jg y^j-y gniall. This implies that the 

Yafi in (69) are very small. So we expect that the Ya^ are also small and that we can set 
them to zero in all the convolution terms of the TEA equations. In this limit, the energy is 
given by inserting the asymptotic form Y^^o, (69), in the expression for the energy (71). One 
would be tempted to expand the logarithm in (71), since Y^^o is very small. However, Ya_o 
has a double pole a m = 0, or g = 0, coming from the boundary dressing phase. In other 
words, it behaves as 

Y.,0 ~ ^ + 0{1) (73) 

q^ 



^ The indices of Yi,m or y„i,i run over m > 2. For Yt^o they run over b > 1. The same as in 
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for small q. We can then write the integrals in (71) as 



oo 



Q \ q^ ) 

In the second term we can certainly expand to first order in Y^o and G^, which produces a 
result which is of order e-2£^m(i+i)_ ^jj-^g gj^g^^ term, however, gives Ga/2 ~ Q-Ern{L+i) ^ which 
is bigger. So we get 

-. oo 

^-"sE^" (75) 

a=l 

But this is precisely the same as what we got in section 2. Namely, (32) is the same as (75) 
after we realize that Ga defined in (73) is essentially the same as (41), using (69). This is 
not too surprising since [23] derived (32) by appealing to TBA equations. In summary, (75) 
agrees precisely with (42). 

In the next section we will perform a weak coupling check of the equations. We will 
derive a simplified set of equations that describe the small angle limit ^, ~ and we will 
expand and solve the resulting equations up to order g^. 



5 The near BPS limit 

When (j) = 9 the Wilson loop is BPS and the energy vanishes. As we deform the angles away 
from this supersymmetric configuration, the energy behaves as 

rcusp(0, e) = -(02 - 9')-^bC\) + oi{<p' - e^f) , a = a(i - ^) . (76) 

The function 5, also known as the "Bremsstrahlung function", is related to a variety of 
physical quantities [8,9]. It was computed exactly in [8,9] using localization. In the planar 
limit we get 

1 71/2 (V^ ) 

4^' /i(\/a) 

On the one hand, this allows us to test the BTBA equation to high loop orders by penetrating 
deep into almost all parts of the equation. On the other hand, the simplicity of (77) suggests 
that, in the near BPS limit, the BTBA equations can be drastically simplified. The equations 
we will find in this limit are not that simple. We hope that understanding how to simplify 
them will teach us how to simplify TBA equation in general. 

In this section we will study the BTBA equations in this limit. We will show that the 
BTBA equations can be reduced to a simplified set of equations. We will then solve them 



g = ^^ '^^ ^ +0(1/Ar2) (77) 
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to 3-loop order. Here we restrict the discussion to 9 = 0^, so that A = A and Tcnsp{4',(^) = 
— 0^-B(A) + (9(0^). We also set L = to extract the cusp anomalous dimension. It is 
important to note that now (p is the smallest parameter. In particular, it is smaller than A. 

In this small angle limit, the momentum carrying Y-functions are of order Ya^ = 0{(f)'^) 
and therefore very small. This limit reminds us of the large L asymptotic limit where the 
momentum carrying Ya^s are exponentially suppressed. However, as opposed to the asymp- 
totic limit, in the small angle limit, we cannot drop the convolutions with the momentum 
carrying Fa.o's- Instead, we remain with a simplified set of non linear equations. The reason 
is that the large value of log Ya^o is not due to the sources in the BTBA equations. Instead, 
it is due to the fact that the fermionic F-functions {Yii and 12,2) approach —1 and lead to 
a big contribution through the log(l + Yi^i) and log(l + ^2,2) terms in the convolutions. 

5.1 The simplified equations at small angles 

As the momentum carrying Y-functions are small, they only contribute to B{X) through 
their double pole. We define C^ as the coefficient of the double pole at u = 0, 



hm Yafi 



2u " 



(7J 



The energy, which is dominated by the value of Ya^ at the double pole, reduces to 

^ = ^ V '^" (79) 

where square root factor comes from the g — )■ limit of {q/2u), see (70). 
In this small 0-limit, the other ^-functions can be expanded as 

1^1,1 = -1-0^^ + 0(0"), y^,i = 3^^[l + 02(l]„-A'^)/2]+O(0^), (80) 

where we assumed that to leading order Fi 1 = Y2,2 = —1 and Fm,i = ^i,m- It is not difficult 
to see that this assumption is consistent with the BTBA equations. Moreover, we find that 
the functions ^2^ drop out of the equations. 



®The general near BPS case, with 6* 7^ 0, has a similar degree of complexity. In fact, we have explicitly 
expanded the equations up to second order in A and verified the corresponding expansion in (77). But we 
will not give the details here. 
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We find that the BTBA equations (62)-(69) reduce to 



^ = - + K^^i * 



$ 



+ K, 



m—l 



t\^rr 



Xr 



y. 



i + ym 
i + ym 



+ 



+ 



x„ 



-Km-l,n-l * log (1 + yn) - K^~l* log 



7^C,7^i°J^(M,0) 



m 



-t^m—l,n~l * 



X -^^ + ^ 



TrCr 



7?(oi) , »(oi) 

'^mn ^^ *-^m— 2,n 



A. 



.(10)' 



^ 



?(io): 



$ 



Kr*iog(^)-srr*iog(^) + 

C = (-l)VF(a,<?)^e^— 



■-r,(lo) , »(10) 



* log 





(81) 




(82) 




(83) 




(84) 


M,0) 


(85) 


.v,.,)l 


(86) 
«=o 



(87) 



±a] 



where z^^""^ denote the values of z^^"^ at g = (25). In (86) we are evaluating the non- 
convoluted variable of the kernels at m = 0. The hat on * is a convolution over the range 
|m| < 2g. F{a, g) is given in (40). These equations are derived by implementing the expansion 
of Y-functions (125) in the TBA system of equations (62)-(66). Let us make a couple of 
comments. First, the factors of 1/2, 1/3, m'^/3 

stand for the subtraction of the asymptotic solutions. These read 

ym = — o 7 5 ^m = 5- (88) 



^ = <l> 



m^ — 1 



Second, note that in the EPS vacuum where = 0, the TBA equations are not well 
defined and need a regularization. A regulator commonly used is a twist for the fermions [41]. 
Here, the angle can be viewed as a physical regulator. As opposed to other regulators, the 
leading order solution 3^^ is a non trivial function of the coupling. 

5.2 Weak coupling expansion of the small TBA 

To test the BTBA equations, we have solved the small angle simplified equations, (81)-(87), 
up to three loops. In this section we will present the results. The derivation is given in 
appendix F. 

The small TBA equations, (81)-(87), are certainly simpler than the general TBA equa- 
tions (62)-(66), but they continue to be non-linear. However, if we make a weak coupling 
expansion we obtain a linear system of integral equations order by order. 
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To solve these linear equations we find it useful to first simplify the TBA equations as 
in [42,43]. To simplify (81) and (82), we take a convolution of the equations with s * s^^ 
where 



s{u) 



(89) 



2 cosh(7r'u) 

The other equations can also be simplified as shown in the appendix F. Then (81)-(87) 
become 



X2 



$ + ^ 



-2s* 



,(01) 



1+3^2 



+ 2ns * n'Z'iu, 0)Cn -nCa Ka{u, 0) 



log ym = S* Im,n log " + Sm,2 ^* log — 



Ar, 



S* L 



x„ 



1 + 3^. 



+ 7r5C™ + (5™,2 5*($-^) 



A. 



>{10)' 



Kr*iog(^)-C°^»iog(^; 



-7?(10) , »(10) 



* log 






C, = (-l)VF(a,^)^e^— 



(90) 
(91) 

(92) 

(93) 

(94) 
(95) 



where /m,,n = 5m+i,n + ^m-i/ra and Kj/^a is defined in appendix F. 

Now expanding the functions \1/, $, 3^„ and Xn in powers of g^ ^ we can obtain them order 
by order by solving a linear system of equations. Up to three loops (see appendix F for 
details) we find that 

' -4 47r2 20' 



a„2 



Ca=4(-1)'^(7 



-\Y 



TT 



^^ + 16(-1)" 



vr 

y 



a2 + a4 



^•^ + 0(<7« 



Finally, the relation (79), we obtain the expression for energy up to 3-loop order^ 



E 



/-^^^+/^ + 0(/) 



A 



A^ 



+ 



A3 



167r2 3847r2 61447r2 



+ C(A^ 



In perfect agreement with the expansion of (77). 



6 Conclusions and discussion 



(96) 



(97) 



In this paper we have considered the problem of computing the quark anti-quark potential 
on the 3-sphere in A/" = 4 super Yang Mills in the planar approximation. Since the planar 



We encounter the sum X)^i(~l)° 



This can be understood by regularizing it as 



hm 



0-i-O 



a=\\ ^) ad) 
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theory is integrable [12], we expected to be able to derive an exact expression. Indeed, we 
found a system of boundary TBA equations (62)-(66) which determines the potential as a 
function of three parameters: the planar coupling A, the geometric angle 0, which sets the 
angular separation on the 3-sphere and an internal angle 9 which is the relative orientation 
of the coupling to the scalar field for the quark and the anti-quark. 

This quark and anti-quark configuration gives rise to an integrable system with a bound- 
ary. This is most clearly seen in the string theory picture where we have a string going 
between the two lines on the boundary. One might be surprised that we have a boundary 
since the string is infinitely long. However, note that the local geometry of the string near 
the boundary is AdS2i which indeed has a boundary. The energy is then the ground state 
energy , or Casimir energy, on the strip and it is given in terms of the solution of the TBA 
equations (71). This is the energy of the flux tube connecting the quark and anti-quark. 
These TBA equations should also enable one to compute the energies of excitations of the 
flux tube. These correspond to operators that are inserted on the Wilson loop. 

The quark anti-quark potential on S^ is the same as the cusp anomalous dimension as a 
function of the angles, Tcusp{,(t), 9, A). 

The derivation of the boundary TBA equations is similar to the one in other integrable 
models with boundary [18]. A crucial step is the determination of the boundary reflection 
matrix. The matrix part is fixed by the symmetries and the dressing phase was found by 
solving the boundary crossing equation and the final answer is in (17), (18). Since there is 
always a certain amount of guesswork in determining the dressing phase, we have checked 
it at strong coupling and we have seen that it gives the right value both in the physical 
and mirror regions. A crucial feature of the dressing phase is that it contains a pole at zero 
mirror momentum. This is crucial for the proposed phase to work at weak coupling. Note 
that the boundary dressing phase is responsible for the leading order contribution in the 
mirror picture, while it only starts contributing at three loops for anomalous dimensions in 
the physical picture. The pole simply means that the boundary is sourcing single particle 
states. 

The BTBA equations were written in (62)- (66). They look very similar to the bulk TBA 
equations [34,36,42], except that the boundary conditions for large u are different. They 
now depend on the angles. In addition, for the momentum carrying nodes, the Ya^, there is 
an extra source term involving the boundary dressing phase. 

We have obtained a simplified set of equations, (90)-(95), which describes the small 
angle region, 0, ^ <^ 1. In this region, the simplest way to solve the problem is through 
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supersymmetric localization, as explained in [8]. The planar answer is 



.20 , .o..4^ o 1 VXhiVX) 



4vr^ /i(VA) 

So, we know the answer by independent means. Thus, these simplified BTBA equations 
should reproduce (98). Indeed, directly expanding these simplified equations up to third 
order in the coupling we reproduced the expansion of (98). However, these "simplified" 
equations are vastly more complex than the simple Bessel functions in (98)!. Thus, there 
should be a way to simplify these equations much further and directly get the simple answer 
(98). Hopefully, the methods used to simplify the equation will also be useful in order to 
simplify the full BTBA equations for general angles. Note that in [44] the TBA system for 
closed strings was reduced to a set of equations involving a finite number of functions. It is 
very likely that the same method works in our case. 

Notice that the simplified small angle equations connect the integrability and the localiza- 
tion exact solutions. In particular, computing the function B by both methods would enable 
us to see whether the coupling constant A that appears in both approaches is the same or 
not. Of course, we expect them to be the same for TV = 4 super Yang Mills. However, if one 
could generalize the discussion in this paper to Wilson loops in ABJM theory [45], then this 
small angle region could enable us to compute the undetermined function h{X) that appears 
in the integrability approach to the ABJM theory [46]. 

In principle, one might wonder whether the Wilson loop leads to an integrable bound- 
ary condition. We have found that the reflection matrix obeys the boundary Yang Baxter 
equation. The TBA equations were derived assuming integrability. So all the checks we 
performed on them are further evidence that the Wilson loop boundary condition is indeed 
integrable. 

There are further checks of the equations that one should be able to do. In particular, 
one would like to reproduce the BES equation [6] for y9 — )• 00. 

It would also be nice to take the small 6 = n — (p limit. In this limit the answer should 
go like 1/6 and probably one can obtain again a simplified equation for the coefficient. This 
determines the quark anti-quark potential in the fiat space limit. 

One should also be able to take the strong coupling limit of the equations and reproduce 
the result derived from classical strings in AdS^ x 5*^ in [47, 48] . It is likely that the ideas 
in [49, 50] would enable this. 

Though solving the TBA equation analytically looks difficult, it should be possible to 
solve the equations numerically. The problem should be very similar to the one solved 
in [51]. 

It would also be nice to study the problem of determining the open string spectrum on 
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the AdS^ X S"^ or AdS2 x S^ D-branes which also preserve the same amount of symmetry. 
The only difference with the current paper should be a different choice for the boundary 
dressing phase. For this reason, the TBA equations would be the same, except for the choice 
of the boundary dressing phase. 

The study of perturbative amplitudes at weak coupling has found remarkably simple 
underlying structures. It would be interesting to study these structures in the context of 
the cusp anomalous dimension, where we have a function of a single angle (p. In particular, 
it would be nice to see how to connect those structures with the TBA approach described 
here. This would most probably lead to both a simplification of this TBA approach as well 
as some hints on the exact structure underlying the amplitude problem. 

Throughout this paper we have considered the locally BPS Wilson loop which contains 
the coupling to the scalar, as in (2). Of course, one can also consider the Wilson loop which 
does not couple to the scalars, W = trPe^^^. It would be interesting to see whether this 
leads to an integrable boundary condition. At strong coupling this loop leads to a Neumann 
boundary condition on the S^ [52], which is classically integrable. ^° 

Note: We were informed that similar ideas were pursued in [24]. 

Acknowledgements 

We would like to thank N. Arkani-Hamed, B. Baso, S. Caron-Huot, N. Drukker, D. 
Gaiotto, N. Gromov, I. Klebanov, P. Vieira and A. Zamolodchikov for discussions. 

A. S. would like to thank Nordita for warm hospitality. This work was supported in part 
by U.S. Department of Energy grant #DE-FG02-90ER40542. Research at the Perimeter 
Institute is supported in part by the Government of Canada through NSERC and by the 
Province of Ontario through MRI. The research of A.S. has been supported in part by the 
Province of Ontario through ERA grant ER 06-02-293. D.C would like to thanks IAS for 
hospitality. The research of D.C has been supported in part by a CONICET-Fulbright 
fellowship and grant PICT 2010-0724. 

A Reflection Matrix 

With the conventions we are using, the (canonical) diagonal symmetry generators are 

-^D + = L\ - L+j^ , R-D f^ = R\ + R\ , Qbh = Q aT iQ\ ) 

L^^ = L%-L\, S^^ = S\±iS%. (99) 

The generators (99) give rise to the diagonal 5^(2 |2)d residual symmetry. Now we should 
determine how a bulk magnon transforms under the diagonal iM(2|2)D. A bulk magnon 



•"See [59] for a systematic study of classically integrable boundary conditions. 
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transforms in a representation {\Zl(a,b,c,d)y 0{a,fe,c,d)) of the bulk symmetry su{2\2)]^ x sm(2|2)^. 
The quantum numbers 

tC ( X \ 77 X ( X^ \ 

a^Vsn, '' = v^f(--ij. -=-v^^. '' = -Vs-(^-i). (100) 

characterize the action of the fermionic generators. For the left fundamental (0^, 0^, ip^,ip^) 

Qll^lj^) = be'^^eM"), Sll^P^) = d6^jr), (101) 

and similarly for the right generators Ql^ and S"*^ and on a right fundamental (0^, 0^, ?/'"'", ^~). 
C is a phase and unitarity requires |?7p = i {x~ — x~^). 

The left part of a bulk magnon, \Zl(a,b,c,d) is also IZiP ?, ^ rf), with identical quantum numbers 
However, the right 0(a,6,c,d) needs to be re- arranged to transform canonically under the action 
of the diagonal symmetry generators. That can be achieved defining, 

(0\ 0^ ^+, ^^) := (0\ 0^ -7^^, +2^+) , (102) 

which turns out to be a IZlP^ _^ -cd)- -^^^ ^° ^^^ signs in —b and — c we should interpret the 
right part as a magnon with quasi-momentum —p and phase (e^^ under su(2\2)-£,. Therefore, 
with the change of basis (102), the original left and right parts of the bulk magnon transforms 
in the following tensor representation of sm(2|2)d 

0(a,6,c,d) ®\Zl(a,-b,-c,d) = V (p, C) ® V (-p, (e'^) . (103) 

Diagonal 5^(2 |2)d is preserved during the reflection, which fixes the boundary scattering 
matrix up to a phase factor. We should take into account that the multiplet labels change 
with the reflection according to, 

V (p, C) ® V (-P, (e'P) ^ V i-p, C) ® V {p, Ce-*^) . (104) 

The reflection matrix TZ-R,{p, — p) intertwines the same representations as a bulk S-matrix 
S{p, —p), and therefore the two must be equal up to a phase [16]. The resulting reflection 
matrix is given by 



7^R 10; X 0i^) = AnipMil X 0j>) + BnipM% X 05I) + lCKip)e'%^\r_^ x 4,^^) , 

nn li^i X V^!p) = Dr,(p)|V^5 X V^f ) + Er(p)|^!!; x ^|1) + iFR(p)e,^6^'^|0ip x 0j) , 

7^R|0jx^!p) = GR,(p)|^^,X0P + i7R(p)|0ipX^f), 

7^R|V^Jx0ip) = KnipM%x4>l) + LM\^-pXi'i)- (105) 
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where 



X+ T]ir]2 



2x+ rji 2x+ 772 

77j are the values of rji after the reflection. The choice for rj in the so-called string theory 
basis is 

1_ vp_ 



■r]{p, C) = C^c 4 \/ix~ — ix^ , (107) 



B Solution to the boundary crossing equation 

In this appendix we solve the crossing equation for the boundary dressing phase 

x~ + -^ 

(^B{p)crB{p) = , , \ , (^b{p)(^b{-p) = 1 • (10^ 

x+ + -V 



We follow a procedure similar to the one described in [21, 22]. First we introduce the 

Zhukovski variable u, 

1 _ u 

x{u) g 
so that the crossing equation can be written as 



x(M) + --y-^ = -, x^:=x(m±|), (109) 



u — - 
as{u)al{u) = ^ . (110) 

u + 



2 



The index 7 means the analytical continuation along a closed contour 7 that crosses both 
x^{u) and x~{u) cuts. We are also going to assume that 

a^[x^[u),x^[u)) = ^^^. (Ill) 

-d 

In terms of the shift operator D := e2 ", the crossing equation reads 

{G{x{u))G{\lx{u))f-'''' =M^"'-^. (112) 
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Therefore, our crossing equation is of the form 

G{x{u))G{llx{u)) = M^(^) . (113) 

Naively one would say that F{D) = —1, but that would associate the cuts of x{u) to 
G{x{u))G{l/x{u)). Instead, we will use 

F.P) = ^3^ + ^3^ = E^" + E^""' (114) 

n=l n=l 

for some integer k. Different values of k would lead to different expressions for u^^'^^\ A 
posteriori we will analyze what values of k are consistent with the crossing condition. With 
this Fk{D), we get 

^FUD) ^ g^p (f^^Ogiu^ + ^) J . (115) 

As in [22], we regulate the divergent sum in (115), by taking the derivative of the exponen- 
tial's argument, performing the sum and integrating back. As a result, up to an irrelevant 
integration constant, we obtain 



f« r(i + ^)r(i-^)- ^'''^ 

We should now check consistency with original crossing equation, so we compute: 

(117) 



smhCfu)\'' ^"_sinh(f(«+|))«-| 



fu ; sinh(f(M-|))M+2 



Consistency with the crossing condition (110) requires 



sinh (f (m + D) _ 

'27r /„, i\\ ~ -^ ' 



1, (118) 



sinh(f(«-|)) 

and we then choose the value k = 1. 

We still need to solve for the function G{x). Let us define G{x) := e^^^^\ The crossing 
condition (113) imposes 

/ / NN / / NN 1, /sinh(27rM)\ 

xixiu + tO)) + xixiu-tO)) = -Aog( ^^^ M . (119) 

The kernel introduced in the bulk case can also be used to solve our problem. Indeed, the 
kernel 

2g+i0 , 

/dvi x[U) r^ 1 
1^ , ^y fiw) , (120) 
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satisfies 

{K ^f){u + tO) + {K^f){u- 20) = f{u) , (121) 

if 1^1 < 2g. Tlius, equation (119) is solved by 

X{x{u)) = -iK ^ log {-^^—1\ . (122) 

Up to a term that cancels when we compute the difference xi^^) ~x{x~)i w^ can write x{x) 
as a contour integral, 

r dz I ^ /sinh[27r^(z+i)]\ , , 

X{x) = $ X =-i i — log _^ (\i^ ■ 123 

J 2'nix-z \ 2Tig[z + ^) J 

\z\=l 

This 'DHM' representation [53] of the solution for the crossing equation is valid for |x| > 1. 
As X moves towards the interior of the unit circle, the function defined by the contour integral 
is discontinuous, as it picks up a residue at x. Since we want the analytic continuation to 
the interior of the disk to be continuous, for |x| < 1, we have instead 

/sinhf27ro(x+ -)1\ , ^ 

X{x) = <l>{x) - 1 log ( 2lg{x+'' 1 ■ ^^^^^ 



X ' 



This analogous to the analysis of [54] for the bulk dressing phase. As x moves inside the 
unit disk, some of the branch cuts of the logarithmic term in (124) can be crossed. There 
are infinitely many of those branch cuts, corresponding to the zeros of smh.{2Tig[x + -)). 
However, crossing such branch cuts could only produce 27r terms in x which are, in any case, 
irrelevant to the boundary dressing factor o"5(x^,x^) = e*'^*-^ )-«x(a; )^ 

Being careful about the crossing contour one can check that the final a^ obtained after this 
procedure indeed solve the crossing equation (108). The unitarity condition is a consequence 
of the fact that x{~x) = x{x) plus the the fact that as is a ratio of a function of x~^ and 
the same function of x~ (111). 

We can expand the contour integral $(x) in negative powers of x for large |x| values. 

The coefficients Cr{g) can be expanded either in the weak or in the strong coupling limit. In 
the weak coupling limit we obtain 

c.to)-E"7""^T"''"^";'- (126) 

^ r(l-2n)r(-|+n+l) ^ ^ 
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Notice that for odd r the coefficient Cr{g) is vanishing, while for even r is order g^ . 
In the strong couphng hmit we find 

cJg) = 2gi / dte^*n costl + Oil) = — ^-—^ — ^ + Oil) . (127) 

J 1-r^ 



To evaluate the contour integral for |x| < 1 we can use the identity 

$(a;) + $(l/x) = $(0) , for|x|^l, (128) 

and the expansion 

u.(2n-l)!! 

n=l 



$(0) = -^Y: 9'''i-lQr^-^^^ ■ (129) 



C Luscher correction at strong coupling 

In this appendix we consider an open string operator of the form BiZ^Bf{6, 0) and compute 
the leading correction to the energy for large L, this correction goes as e-(constant)L^ -^^ ^jjj 
compute the correction for = 0, ^ 7^ at leading order in the strong coupling expansion. 
In this case we have a string that moves on AdS2 x S^. It is convenient to decouple the 
AdS2 and S^ problems by choosing a worldsheet gauge where —T^^^ = T^_^ = 1. We fix the 
solution on the AdS2 part. This AdS2 solution is completely characterized by the extent of 
the spatial worldsheet coordinate a, which we take to run between [— s/2, s/2]. The other 
worldsheet coordinate is r. In particular the spacetime energy A of the solution is fixed, 
once s is fixed. As we vary the parameters we will see that L will change, 6 will change, and 
so will A — L. 

So we now concentrate on the solution on the S^, which we parametrize as 



Xi 



+ %X2 = e'^Vl-p'(^) , 2:3 + tX4 = p(a)e*'^('^) (130) 



Inserting this in the Euler Lagrange equations for the string and imposing the Virasoro 
constraints, T±± = 1 one finds two integrals of motion, i and 7. They are given by 

i = /)V (131) 

and 

(^ = -f - (7^ - 1)P' + I'p' (132) 

1 — p^ 

The boundary conditions are p'(0) = 0, p(s/2) = 1. Let us define po to be the value of p at 

a = where the derivative vanishes. It is a root of 

= -f- (7^ - 1)P^ + I'pt (133) 

38 



By using (132) we can write the following expressions 

i = r^^ /rA/n (134) 

9 f^ i 

2 = / ^^ /^ 2/n (135) 

- = 2g da-f\xi + ix2\'^ = 2g'y dp y= — (136) 

'^ Jo Jpo VIJ 

D = -^^-(f-l)/ + 7V = (p^-p^)[7^(p^ + P^)-(7^-l)] (137) 

From the first two equations we should find po and 7 as a function of s and 6, and then we 
can find the expression for L and for the energy. We want to find a solution where L is very 
large. 

This happens when po ^ and 7 — ;■ 1 and £ — )■ 0. More precisely, we need to scale them 

as 

7 = l + e/2, £ = e^, p = v^t; (138) 

where f is a new rescaled variable and i is fixed as e — )■ 0. Now, to leading order in e we find 
that (133) becomes 

= -^-^o' + ^o, or vl = \^ (139) 

The integral for ^, (135), becomes negligibly small away from p ^ p^ since there a factor of 
i multiplying. So it receives all its contribution from the small p region, namely the finite v 
region, see (138). We can write 



Q 



tan^ (140) 



2 2 j^j, ^^^ 

b = {v^-vl){v' + vl-l) (141) 



We can similarly compute the integral for s. 



1^ 0'vh^i7B-?)^0^wh <"^' 

dvvl^-\]+\og2- logpo (143) 

\^/D ^ y 

'- = log4-Jlog(^), ^^ = l6e-^ (144) 

2 2 cost^ cost' 

where we used po = y/evo and the result (140), and the definition of D in (141). Here we 
have split the integral in two terms, the first receives contributions only form the small p 
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region and the second, which can be done exphcitly with no need to take the small po limit 
(though we quoted here only the small po answer). 

We now want to compute L. We will compute instead 
L s 



4g 2 



L _s 

4g~2 




(145) 
(147) 



(148) 
l + cos^4e"' (149) 



L s _ 

Ag~2 

L-2gs = -4:g + 16gcosee~^~^ (150) 



Here we have split the integrals having in mind that we want an accuracy of order e. The 
first has an e in front and we made sure that only the small p region contributes. In the 
second we made sure that only the small p region contributes up to order e. The last can 
can be evaluated exactly and we quoted here the small po result. 

Here we have in mind keeping s fixed as we change 6. Under these circumstances A stays 
fixed, since the ^4^5*2 part of the solution would always be the same. In addition, we know 
that for ^ = the result should vanish due to the BPS condition. Thus we find that 

16 L 

A-L = g(l-cose)—e'^9 (151) 

If we changed the angle in the AdS part, then instead of 1 in (151) we would get some 
function of 0. However, since we know that for 6' = we should get zero due to the BPS 
condition, we conclude that for generic angles we get 

16 _A 

A - L = 5f(cos0- cos6')— e "^ (152) 

C.l Strong coupling expansion of the function F 

In order to compare this to the expected answer from the Luscher type correction we need 
to evaluate the function F in (40) at strong coupling. This involves evaluating the function 
$ in (17), (19) at ^[=^"1 at g = 0. When g = we have that 

- l/^[-'^l(0) = ^[+'^1(0) = t (/ITWIieTy- a/(4^)) =tfl + ^ + ...\ (153) 
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which is very close to i, where the strong couphng expansion is tricky, since we cannot use 
(22). We need to compute 



logF = i^{y) -i^{l/y) 



vr 



dt- 



(y' - 1) 



[1 + y^y - 4y^ sinH 



log 



sinh Ang sin t 
Ang sin t 



(154) 



with y = 2;["1(0). Then the y dependent factor can be well approximated by 



iy' - 1) 



=xW(0) 



4^ 



smH + 4^, 



(1 + y'^y - Ay'^sin^t 
We now insert this into the integral (154), and split the integral into two pieces 
log F = ri + r2 

1 



n 



r2 



r2 



27ig 



27ig 



dt 



dt 



2+ I _ai_ 
16g2J 



sin t + 
1 



2+ , _g?_ 
I652J 



sin^ t + 



Ang sin t = 2a log 
log 



a 



+ o{l/g) 



(155) 

(156) 

(157) 

(158) 



Aa 
a^ 2 I ^ 2 2 ^°S 



1 - e-" 



Sirg sin t 
= 2[aloga-a-logr(a + l)] (159) 



where we have defined t = v/{87Tg) in the integral for r2 and taken the (7 — ?■ 00 limit. 
Summarizing, we get that the leading strong coupling approximation is 

06a „2a 



F{a,g) 



e^Hal)^ 



(160) 



D Evaluating t{q) 



Let us now evaluate t{q), given in (31). In order to perform the trace over the matrix indices 
it is convenient to write the refiection matrix in terms of the bulk S matrix. The reason is 
that we will be able to use the bulk crossing equation to simplify the form of t{q). We start 
by writing 



K 



1 + 



(2+)2 z+ + l/z- 



^^'^^{q) = a,{z\z~)—^J \ 7 7_ g^g(-p,p)TiT#C^^C^^ 



1 + 



{z-y 



z- + 1/z- 



(161) 



where S is the full bulk 5* matrix for one of the sm(2|2) factors, which obeys the crossing 
equation, with the identity in the right hand side. Here p is the value of the momentum 
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analytically continued so that we have x^ in the mirror region. The square root arises 
because of a conventional way to define the bulk S matrix. It is cancelling a square root 
in a the phase factor of the bulk S matrix. The matrix T converts the dotted indices into 
undotted indices. It arises in the precise implementation of the "unfolding" trick, where we 
replace a bulk magnon with momentum p that transforms under iM(2|2)^ into two magnons 
of su(2|2)£i, one with momentum p and the other with momentum —p, see appendix A. We 
then see that K in (30) is essentially the same as K, but evaluated at —p. More explicitly, 
we can write 

^aaM^) = -. (-^, -^) [${5|^^^$j^«(l>, -m^T-CoACHN . (162) 

We can now insert this into the expression for t{q). We will need to use that 

^^B(-P,P)S^C^^C^^5ff,(p, -p)J:'JCgaChn = Tr[S]2 , (163) 

here E^ = T^T^, where 

S = diag(l,l,-1,-1) = (-1)^. (164) 

This arises because the action of the charge conjugation changes in the basis given by (102). 
The charge conjugation matrix can be taken to be 

Cs = ( -"; ° ) . (165) 

Equation (163) arises from the repeated use of the crossing equation. Ignoring charge con- 
jugation matrices the identity we need is 

S{-p, p)S{p, -p) = S{-p, p)S{-p, p)S{p, -p)S{p, -p) = 1 (166) 

where S denotes the full S matrix. When (163) is used, we get (33) in the main text. The 
series of operations we have done are most clearly summarized by the figure 9. First we do 
the unfolding trick. Then the use of the crossing relation (166) amounts to moving the lines, 
as in figure 5, and untangling them. If we introduce the rotation matrix m, we can do all 
the same steps but we have insertions of the matrix m or m~^ along some lines. This is 
represented in figure 9 by the solid circles. Once we untangle the lines as in figure 5, we get 
get an insertion of m on one line and an insertion of m^^ on the other, leading to (34) (since 
the trace of m or m'^ are the same). 

This gave us t{q) for the fundamental mirror magnon. For the bound states we can use 
the fusion procedure. Given the form of a^ in (17), which involves a ratio of a function of 
x~^ and x~ , then it is clear that the fusion procedure gives a as which is the same ratio, but 
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o ■ 

o : 



(c) 

Figure 9: (a) Original picture. White circles represent rotation matrices, one is m and 
the other is m^^. Solid circles represent charge conjugations, (b) Unfolded picture. In the 
unfolding of the dotted indices there is a change of basis that produces the S's. (c) Untangled 
picture . After using crossing we get two independent traces of the matrix m. 

evaluated at 2'^"' and 2;[^"l The corresponding matrices are determined also by multiplying 
the matrices of the elementary constituents. All the manipulations we used above can be 
used again for these matrices. In particular, we can untangle the lines as in 9. The only new 
thing we need to understand is the set of states of a magnon bound state and the action of 
the matrix m. The matrix m is 



rua 



diag(e 



ie-i{a-l)<l) _ _ _ ie+i{a~l)<l) -ie-i{a-l)<j> 
a+1 a— 1 



-ie+i{a-l) 



^ia<i> g-i(a^2)0 



J{a~2)4>\ 



(167) 



This can be understood as follows. The mirror magnon bound state arises from a SL{2) 
sector fundamental magnons, giving rise to an SL{2) representation of spin 2j = a, these 
lead to the components with the a + 1 bracket. The other elements arise from acting with the 
supercharges in sm(2|2). Then we see that the trace gives the result quoted in (36). Thus, 
the cancelation of the bulk S matrices, plus the form for as for the bound states, together 
with (167) lead to ta{q) in (35). 

E Derivation of the BTBA equations 

In this appendix we derive the BTBA equations presented in the main text (62)- (69). We 
will do so in two different ways. In the first way, presented in section E.l, we will follow 
the original derivation of the TBA for the spectrum [37, 38] by embedding our open ABA 
equations into the closed ones. In the second way, presented in section E.2, we will take a 
more direct route and derive the BTBA from the thermodynamics in the mirror picture [17]. 
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E.l Derivation by embedding into the bulk system 

The spectrum of closed strings or single trace operators in A/" = 4 SYM is governed by the 
so called Y-system [37]. The Y-system is a set of functions Ya^s{u) characterizing the ratio of 
the density of bulk excitation to the density of holes in the mirror channel [31]. The indices 
(a, s) stands for rectangular representations of the bulk S't/(2,2|4) excitations and u is the 
spectral parameter. The Y-function are subject to the general functional equation 



ya+l,sYa-l,s (1 + ^a+l,s)(l + ^a- 



1,^; 



where /^ = f{u± if^)- For the bulk excitations, the Y's live in a fat hook bounded by 
Fo.s = oo, 1^2,|s|>2 = oo and Ya>2,±2 = 0, (see fig 8. a). The Y-system is also equivalent to the 
Hirota equation for the T-functions as 

T T 
Ta,sTa,s = Ta,s+iTa,s-i + Ta+i,sTa-i,s , where K,, = "'"^^ °''~^ (169) 

-^ a+ljS-' a—l,s 

and enjoy a gauge invariance under Ta^s — ^ 5''^'^^'^^^a,s- For more details, see [37,38]. The 
TBA equations for the Y functions are the solution to the functional relation (168), subject 
to the relevant boundary conditions and analytic behavior. The procedure of deriving the 
TBA equations in this fashion was carried out in [37,38]. 

Operators on Wilson loops are associated to a string, or a spin chain, with boundaries. 
Suppose we start with such operators and go to the mirror picture where space and time 
are interchanged (see figure 2). In the mirror picture, one have exactly the same system of 
mirror particles as in the closed case. In the limit of large T, see figure 2, these particles 
live on a large closed chain governed by the mirror asymptotic Bethe equations [35]. The 
boundaries in the original physical picture are mapped to two boundary states in the mirror 
past and future. These boundary states are determined by the boundary reflection matrix 
as we discussed in section 2. As opposed to the closed case where all mirror states are traced 
over in the partition function, in the overlap between boundary states only a subset of mirror 
excitations are summed over. Moreover, the weights of these excitations in the summation 
lead to a new asymptotic behavior for their densities. We therefore expect the the Y-system 
and the TBA equations to be identical to the ones in the closed case modulo projections and 
new sources. 

In this section we will exploit that relation to derive the boundary TBA equations. That 
is, we will first map the open sm(2|2)d ABA equations (56)-(58) into a folded version of 
the standard closed sm(2|2)^ ones. That map is nothing but the embedding of the diagonal 
sm(2|2)d excitations in the full su{2\2)i^ x su{2\2)-n by restricting to singlet excitations of 
the diagonal sm(2|2)d preserved by the boundary. Having done so, the corresponding BTBA 
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equations will follow from the derivation of the closed TBA ones [38]. We will only add 
the angles that enter as (diagonal) twists and the boundary dressing phase that enters as a 
momentum dependent chemical potential for the momentum carrying excitations. As most 
of the details are the same as in the closed case, we will be brief. 

As we discussed in section 2, and in appendix A, the diagonal su{2\2)d magnon excitations 
transforms in a tensor representation of the bulk sm(2|2)l x sm(2|2)r. Correspondingly, the 
ABA equations can be embedded into the closed ones together with the addition of the 
reflection matrix. That is, (56)- (58) can also be written as in a redundant way as 



1 (g) e^<^-^ 
4 A Ro\u) 



g(-)Q+ 



«i,fc 



B(+)g2 

Qi Q2 Q3 



Qi Q2 Qs 



U2,k 



!'->Qj 



R(+)Q; 



<s> 



6 



18) 



-i(j>+ie 



X 
X' 



"3,fc 

S{u 



R+Bt 



X 



M(+)Q2 
-2ie Qi Q2 Qz 



"3,fe 



-i(j>+ie 



Vl V2 ^3 



-^Qt 



-U2,k 



-«l,fc 



Siu 



RtBt 



"4,fe 



(170) 



where 



Ri{u) 
Riiu) 



K±), 



u 



S{u) 



n 



X[U) — Xij 



\ {^i.y 



2 ' 



x{u) + Xij 



K, 1 



BAu) 



Biiu) 



Ki 

n {x,,y/' 

-^ (x(m) - xj){x{u) + xf) 

i=i 
K4 



K, 



(x^Xj )^/2 



- - xiu) -^'J 
TT x(u) + -^^J 

]n>(±)|',,\ _ 


K4 

11 

i^4 






(s''^7)^/' 



s{u) = l[sMu),-Pjf 



(171) 
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where we used that Sq{pj, —p{u)) = So{p{u), —pj)- 

Note that equations for nodes 1-3 in (170) are equivalent to the ones in nodes 5-7.^^ 
However, here we will think about these as describing the left and right parts of the bulk 
magnon excitations correspondingly.^^ We see that momentum carrying excitations come in 
pairs with opposite momenta. The nested level excitations also comes in pairs of rapidities 
Ms = —Ms, Uq = —U2 and u-j = —Ui. That is, Qi{u) is nothing but Qs-iiu) with roots 

Uij = -Us-lj. 

We will now use that embedding to derive from them the asymptotic values of the left 
and right T-functions. Of course, these are not independent and we only have one set of 
independent ?u(2|2)d T-functions. 

In the asymptotic L — )■ oo limit the right (T^ ,,>o) and left (T^ ,,<o) T systems decouples^'^. 
The asymptotic values of Ti i and Ti _i, whose analyticity would lead to equations in nodes 
1-3 and 5-7 in (170) correspondingly, are derived as in the periodic case [37]. In what 
follows, we will use bold face (Y, T) for the asymptotic values of Y's and T's. Up to a gauge 
transformation, the asymptotic Tf^ and T^_^ are 



R 



rpn 

-^1,1 



rj.L 



r(+) 



Q2Q3 

^ie Q2 Q3 _ 
Q2Q3 



— e 



l-^-)Qt 



JO 






+ e 



Q2 Qi 
"" Q2Qt 

-ie Q2 Qi 

Q2Qt 



o'i-4' 



l+M 



Qi 



-i(f> 






(172) 
(173) 



The right 1-3 in (170) ABA equations are obtained by demanding analyticity of Tf^ when u 
goes to Ui^k — |, U2^k and u^^k + |. The equivalent right equations 5-7 in (170) are obtained 
from the analyticity of T^_^ when u goes to —ui^k — |, — ^2,fc and —U3± + \. In the gauge 



where Tq^ 



-,L/R 
-a.O 



1, the other T functions of the right and left (decoupled) 5'[/(2|2) 



wings are obtained from the generating functional [37, 55] 
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_,,Q2-Qt^-^^^ 



1-e 
-, -1 



id 



Q2 Qi 



-(+) 
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(174) 



as 



>v = ETt 



^D' 



s=Q 



w- 



E(-i)"Ti'r'^"' ^=e" 



idu 



a=0 



;i75) 



^^Remember that in physical kinematics, where (170) is written, jx^"! | > 1 and therefore x^ (—u) — —x^iu), 
x{—u) — —x{u), p{—u) — ~p{u) and e(— m) = e(u). 

^^After rearrangement of the quantum numbers discussed in the appendix. 

^•^Throughout this section the word "asymptotic", means this large L limit, which should not be confused 
with the large u limit. 
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Similarly, for negative s we reverse the sign of the angles and use Q instead of Q. Their 
analyticity leads to the Bethe equations for the bound states. In the vacuum, the T's are 
independent of u. Their constant values read 



,R „T . _x„„, , „, smai 



T^i = T^,_i = (-ir2(cos</>-cose)— -f , (176) 

sm (p 

Sin '■^ f) 
Tf, = T^_,= -2(cos0-cose)^^. (177) 

Using the Hirota equation and the definition of the Y-functions in the right and left decoupled 
wings (169) we read the corresponding values of the asymptotic Y's 



cos^ sin[(s + l)e]sin[(s-l)^] 







Yi,s 


1,1 


COS0 




COS0 


Y,,i 


2,2 — 


n ' 



sin^^ 
sin^ (h 



(178) 



cos6' ' ' sin[(a + 1)0] sin[(a — 1) 

For the asymptotic Yafl (Yafi) we have 



2L 



Ya,o=(^) !yT;TTtiTj, (179) 



Here, ( |p:^ j 4q:^ is a zero mode of the discrete Laplace equation _^ '2" — = 1 [37] . 
It comes about because T^ and T'" in (172)- (175) are written in different gauges. It is 
determined by demanding that li,o(w4,fc) = —1 gives the 4**^ node Bethe equation. We find 

The boundary crossing equation (16) is then obtained by demanding that Yi i is invariant 
under crossing (provided that the bulk dressing phase obeys the bulk crossing equation). In 
particular, in the vacuum we get 

Y,o = 4 , , ,, , ,,^u {^] (cos0-cos0)2?^^^. (181) 

Finally, note that T]^_^(n) can be obtained from T]^^(u) in two steps. First reflect the 
sign of u by considering Ti i(— m). Second, flip the signs of all the shifts and the angles. 
In mirror kinematics, that second step that basically amount to a complex conjugation, 
is equivalent to a gauge transformation [56]. We therefore find that in mirror kinematics 
: T™''(-m) and as a result YJJ;L',(m) = Y'^^{u)}'^ As the closed TBA equation 



^''To sec that Y™Q is a symmetric function note that in mirror kinematics, x'^°'\—u) = —l/x^ °^ and that 
Xi-x) = x{x). 

47 



respects that symmetry, we expect all solutions to the BTBA equations to respect it. That 
is, we have a folded version of the bulk Y-system where Ya^-s{u) is identified with Ya^s{—u) 
as presented in the main text (61). Such a folding is also expected a priory as the Wilson 
lines break the bulk SU{2, 2|4) symmetry (at ^ = = 0) down to OSp{4*\4). For the ground 
state, we expect all functions to be symmetric. 

We conclude that the BTBA equations, once divided by the asymptotic ones, are the 
same as in the closed case. The only two differences are the modified asymptotic solutions 
(178), (181) and the folding Fa,-s(^^) = Ya^s{-u). These are the BTBA equations (62)-(69) 
presented in the main text. 

E.2 Direct derivation of the boundary TBA 

The derivation of boundary TBA equations from the boundary state (27) is fairly standard 
in relativistic theories [17, 18]. Here we need to follow the same steps. 

After doing the flip between space and time, see figure 2, we have a past and a future 
boundary characterized by the matrices K and K which give the probability amplitude for 
creating or annihilating a pair of particles (27), (29). This pair has opposite momenta and 
it is in the singlet representation of an S'f/(2|2)o, since this is a symmetry preserved by the 
boundary state. By independent creation events we can create a multiparticle state. We 
can graphically represent the quantity we want to evaluate as in figure 10(a). Along the 
spatial direction we have a closed circle of length T, we need to solve the Asymptotic Bethe 
equations on this circle and find the Bethe eigenstates. Only the subset of momenta that 
corresponds to Bethe eigenstates can propagate. These states propagate for Euclidean time 
L. These asymptotic Bethe equations are the ones in the mirror theory and were written 
in [32], following [20]. These equations involve various roots, u\,---ui, where u^ are the 
momentum carrying roots. The only new feature is that we are considering states which are 
composed of pairs of particles with opposite momenta. This imposes the condition that the 
^4 roots should appear in pairs. In other words a root M4 should appear together with a root 
— M4. The boundary state is invariant under a diagonal SU{2\2)£). The condition that we 
only have singlets under SU{2\2)i:, implies that if roots «„, with a = 1,2,3 appear, then so 
should roots Wg-a = —Ua- This will be more clearly seen below. 

For this multiparticle state we can perform an unfolding and untangling trick identical to 
the one we did for a single particle in section 2, see figure 10(b),(c). Now we have particles 
for a single sit(2|2) chain, but with insertions associated to the rotation matrix m as well as 
the dressing phase, crsiQ)- We can view these as "chemical" potentials in a thermodynamic 
computation. The Bethe equations in this new picture look slightly more complicated because 
taking a single full particle around the circle, as in figure 11(a), amounts to taking a pair of 
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Figure 10: (a) Original computation. The boundary states create and annihilate pairs 
of particles. The blue and red lines represent the two iu(2|2) representations of each bulk 
magnon. The doted line represents the projection on to states that obey the Bethe equations 
for a chain that is closed along the horizontal direction. The red dots are the rotation matrices 
m in (24) and the black dots are the boundary dressing phase, (b) The same in the unfolded 
picture. We have a single sm(2|2) group. We continue to have insertions of the matrix m and 
the dressing phase. The Bethe equations now involve some operation which also involves the 
bottom part, (c) The same in the untangled picture. Here for each line we take the trace 
over all four states of the magnon. 

particles around the circle, as in figure ll(b),(c). We should project onto states that obey 
these equations. These equations are for the full matrices. We can do the nesting procedure 
and follow the unfolding for the various levels of the nesting. The Bethe equations for the 
various nesting levels can be graphically represented as in figure 12. The final result is that 
the equations are identical to the ones we would get in a situation where we have the full 
mirror theory, but we impose that each root Ua is accompanied by a root —Us^a- 

The final Bethe equations in the mirror picture are 

1 = e'P^'^i^^]JS,,iuiul)S,M,-u'l)llS,3iui,ui)S,si-uiui) (182) 

k I 

1 = llSs,{ulul)Su{ul-ul)l[S32{ulu^) (183) 

k m 

-1 = n^32«, 4)11 ^33(4, «3) (184) 



where Sab are the S matrices (really just phases) between the impurities of the various kinds. 
We have removed the Ui roots by turning them into u^ roots, just to make the equations 
more compact. The Sab are the same as the one that appear in the ordinary mirror Bethe 
equations in [32], see also [20]. In this form, these equations follow in a straightforward 
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Figure 11: The Bethe equations in the various pictures come from demanding that each full 
particle can be taken around the chain producing a total phase of 1. (a) The Bethe equation 
in the original picture, (b) The Bethe equation in the unfolded picture comes from taking 
both a top particle with momentum q and a bottom particle with momentum —q around the 
chain, (c) The same but in the untangled picture. 

fashion from figure 12(c)(d). Note, in particular, that the last factor in (182) can also be 
written as S'43(— 'U4,'U3) = 5'45('U4,M5) with Wg = — M3. So that we can view the whole set 
of equations, (182)-(184), as arising from the full chain, with both iM(2|2) factors, but in a 
configuration with the roots related by ug-a = — ^a, as we discussed above. 

We can simply add the rotation matrices to this picture. They act diagonally on the 
various impurities of the various levels. They appear as chemical potentials. The potentials 
are as follows. The momentum carrying modes, labeled by U4, are in the SL{2) sector and 
have spacetime charges, so we get a factor of e~'^^'^, one e"**^ for the particle of momentum q 
and one for the one with momentum —q. The next root, U3 or Ui, is the action of a fermionic 
generator that changes a fermion into a boson, so that we need to add e*^"'"*'^. Then a bosonic 
generator with e~^*^ is associated to the root U2- Taking into account the number of roots 
of each type in a magnon bound state, we can get the chemical potential coming from the 
angles for each magnon bound state [32]. In addition, the boundary dressing phase also 
appears as a kind of chemical potential. For each pair of momenta q and —q we have a factor 
of 

aBaB=(^B{z^,z')crBi T' +j (185) 

The derivation of the TBA equations for this system follows the same essential steps as 
the derivation for the full TBA for the closed string. There are only minor differences due 
to the constraints on the position of the Bethe roots. For example, for the roots U4 we only 
need to know the density in the region M4 > 0, since there is always another associated root 
at — M4. Thus, all integrals over U4 will run over U4 > 0. We can define the Y functions. 
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Figure 12: (a) Bethe equations for the momentuni carrying nodes in the original picture. 
The doted hues represent the various nested levels, with red and blue for each of the su{2\2) 
factors. Here y represents the rapidity of levels Ui, U2 or M3, and —y is the rapidity of the 
MsjMg or M7 levels, (b) Bethe equation for the nesting levels, (c) Equations for the momentum 
carrying nodes in the unfolded and untangled picture. Note that both the particle with 
momentum q and the one with momentum —q cross the nesting doted line, (d) Bethe 
equation for the nested line. 



as usual, in terms of densities of particles over densities of holes. We then run through the 
usual argument [33], but keeping track of the range of integration, etc. We end up with the 
boundary TBA equations 



log yi,i = te + ict> + Km-i * log ^^^ + 7^S°;) * iog(i + r,,o) 



log 1^2,2 = iO + i(f) + Km-i * log 



1 + Ym,l 

l + Yi^ 

1 + YjnA 



+ Bfl^ * iog(i + r„,o) 



logFi,, = 2i{s - 1)9 - Ks-i,t-i * log(l + ri,i) - K,_, * log 
logy„,i = t2{a - 1)0 - Ka-i,b-i * log(l + n,i) - Ka-i * log 



1 + ^1,1 

1 + ^2,2 

1 + ^1.1 

1 + F2,2 



(186) 
(187) 
(188) 



+ 



-7? (01) , »(oi) 

'^ab "T "a-2.6 



logr„,o = -^2a0 + loglaBCTB] - 2LE^{u) + 



+ 



25„,-7^i") + i3, 



?(ii) 



* log(l + Yi 


),o) 


(189) 


(190) 


+ n,o) + 2 


-r,(10) , »(10) 
'^ab + "a,b-2 


Sy*mMl+n,l) + 



)(10) 



+2Kr * log(l + n,0 - 2Bl\^> * log(l + ¥,, 



(10) 



sym 



sym 



(191) 



^{here) 



where the kernels are listed in appendix E.3 and are the same as in [37,38] but with Y^ ^ 
l/Ya^s^^'^ ■ The values of the angle dependent form of the chemical potential for each Ya^s 
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function follows from the particle content of each mirror bound state associated to each y^ ., 
function [32]. If we set to zero the Y^^a that appear in the convolutions, we get the large L 
solution quoted in (67), (69). These equations look similar to the ones in [57,58] for twisted 
boundary conditions, apart from the folding symmetry and the boundary dressing function 
ctb- 

Finally, the expression for the energy is the one quoted in (71). 

E.3 Kernels conventions 

For the kernels we use the same definitions as in [37, 38] , 



n — 1 7n— 1 



Kn{u,v) = , _ .^ , Kn,rn{u,v)= ^ ^ K2j+2k+2{u,v) , (192) 



rfi + A{u — vY 



71 — 1 Tn — 1 

2 2 



A— "-1 h— m-1 
J- 2 '^- 2 



1 d r{u + iajl + zj, v — ih/2 + ik) 



^-^ ^-^ l-Kidv r{u — ta/2 + ij,v + to/2 + tk) 



Ti — 1 Tn — 1 

2 2 



1 (i h{u + ia/2 + ij, t; - ih/2 + zA;) 



■?= 2^ '^^ 2^ 
5„^(u,i;) = ^Aloga(a;[±"l(M),x[±"^](t;)) , (195) 

r{u,v) = ^ ' : ' , h{u,v) = ' ^ ' ^ ' (196) 

\/x{v) \/x{v) 

Fourier transformations of Kniu, 0) and Kn,m{u, 0) are used at different stages, 

Kniw) = sign(n)e-l"'"l/2 , km,nM = coth M Jg-^l— "I - e-^(™+")] - 5^,„ . (197) 

F Perturbative solution of the small 6 TBA 



Before we start with the perturbative expansion of the small (j) TBA system, let us see in 
more detail how the simplified system (90)-(95) arises from (81)-(87). In order to simplify 
(81), (82) we have to convolute their Km-i terms with s * s~^ and use the identity 

5"^ * K^^l = Ki^rn-l + ^2,m , for s{u) = — r . (198) 

2 cosh(7r-u) 
We then use the TBA equation for X2 and the exact relations 

nf^\u,v) + B^^^\u,v) = Kaiu,v) , (199) 

n^'a\u,v) - B^^^\u,v) = ky,a{u,v) = K{u,v- zf) - K{u,v + 1%) , (200) 
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where 



K(u,v) = ^xl fl ""l ^ . (201) 



Regarding the other TBA equations, we can simphfy them by convoluting with 

The simphfied TBA equations (90)-(95) are useful to solve for the Y-functions in the 
small 0-limit perturbatively. We should begin by expanding the Y-functions in powers of g"^, 

15 

^ = ^(0) + ^(1)^2 ^_ ^{2)^4 _^ 

$ = $(0) + $(1)^2 _,_ ^(2)^4 _^ 

ym = 3^i°^(l + 3^iV + 3^iV + ---) (203) 

The leading orders in the expansion (203) can be obtained by setting setting 9 = Q and 
L = and taking the small limit in the asymptotic solution (67)- (69). Thus we have 

^(o) = $(o)^l^ 3;(0) = _L_, ^W = -^, Ci^)=4(-1)". (204) 

We can obtain the higher order terms by solving the system as follows. By inspecting the 
equations (92) and their perturbative expansions, one realizes that, to any order, the system 
can be solved in this schematic way: ^('^^ - $(*■■) -^ y^^) -^ XJn^ -^ ^(^) + <l>(^) -^ Ci^\ 
provided the order k — 1 functions are known. In the final step of computing Ci (95) 
two kinds of contributions must be distinguished. On the one hand, there is Aconv which 
is originated in from the convolutions with the Y-functions in (94). On the other hand we 
have contributions from the explicit functions in (87) which can be expanded to any order 
independently of the Y-system solution. These give 



a' 



-\ 2 



^0 



[+«] 




a + \l a^ 
16/ 



2 w , .y , , 4 ISOvr^ 630\ 16(7^ , ^^, g. 



where we have used that z^ are given by (25) at g = 0. We kept the first 3 loop orders 
only. 



^^The different choice in tfie ym. expansion is to have similar recurrent equations to those of X„ 
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F.l Small (j) solution at 2-loops 

Let us now solve for the small Y-functions to the next to leading order in the coupling 
expansion. We refer to this as the 2-loop order computation because this order gives rise to 
a correction 0{g'^) to the energy £. Note however, that the Y-functions corresponding to 
this order are 0{g^). 

We start with (90). Using the leading order of C^ and the expression for Ky^a, 

$ - ^ ^ V ^ / ; \ ;. . 206 

^-^ (a2 + 4^2) 

Given that fermionic Y-functions are defined in the interval {—2g,2g), this difference is 
essentially 0{g^). Thus, it will be convenient to use the variable u = u/{2g) that runs 
between —1 and 1. We then obtain 

"^ / 1 '\a o 2 3 

^2gu) - ^(2gu) = 32g^Vl^^Y ^^ + C(/) = ~^^^VT^^ + 0(g'') . (207) 

^-^ a^ 3 

a=l 

This indicates that fermionic convolutions do not contribute to Xn or 3^n Then, we have 

3^«-.*W3^«^^^ = 0. (208) 

The same recurrence equations, thought with different rhs's, will repeatedly appear in the 
weak coupling expansion of the small angle limit TBA equations. In F.3 we present the re- 
solvent of the corresponding recurrence operator (demanding 3^| = and that yL remains 
bounded as m — )■ cxd). In this particular case, we solve (208) with ym = 0. 

The recurrence equation that Xm satisfies is non- homogeneous. It has a non- vanishing 
rhs because of the term its Cm in (93), 

XW - s * I^^^XW^^^I^ = 47r(-l)™. . (209) 



In Fourier space this becomes 

2 cosh p") - l^e, - ^;^«i. - M-ir . (210) 

which can be solved using the resolvent of F.3. Here we write down only the component we 
need 



■2 ~ 

Now, from (91), we have 



^ 3 



2cosh — logf l + e"~j - 1 -e"~ . (211 



$(1) + ^(1) = _3_^ ^ ;t'a) ^ 87r(-l)"5 * K„_i - 47r(-l)"ir„ , (212) 
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where we have used 

7^g^)(«,0) = K,_,iu) + ^K,iu) + Oig') . (213) 

To compute the convolutions we go to Fourier space and get 

$(1) + ^(1) = 4^ log^i + e-kl/2) ^ (214) 

and then 

$(1) + vl/d) = ^ f 1 - -4^) ■ (215) 

m2 \^ sinh(27rM)y ^ ^ 

To compute convolution with this fermionic Y-functions ($ and \E') we shall use the afore- 
mentioned u, at the expense of introducing a (^-dependence which has to be expanded as 

<^'^'\2gu) + vI/(i)(2^S) = ^ - ^/ti^ + 0{g') . (216) 

3 45 

For the next to leading order it is enough to keep only the constant term. Note, however, 

that the second term will also contribute to the next to next to leading order. 

We now compute 



Aeon. = 7^i\°)* log 2^ - i3i\°)* log 2$ + f7^i-) + i3i-i J * log -i±^ . (217) 



1 + 



b2-l 



and we find to the leading order 



271^02 

Aconv = -Y- + ^(9^) ■ (218) 



We are now in the position to evaluate C^ at the 2-loop order. Inserting (218) and (205) 
into (95) we find 

'\' + 0{g') , (219) 



C, = 4(-l)V + 8(-l) 



2 4 



F.2 Small solution at 3-loops 

Convolutions with fermionic Y-functions ($ and \E') start contributing to the next to next 
to leading order. We have already obtained the difference of fermionic Y-functions in (207). 
For the convolution needed in (92) we get 

., ^ 167rV . x^ / o 167rV / x / x 

5* log — ^ s[u — 2gu)*vl — u^ c:^ s[u) . (220) 

Sir O O 

Thus, from (92) using 3^^ ~ ^7?=:t(1 + g^ym + . . . ) we obtain 

3^i^) - s * /.,.3^(^)^^ = ~5^,s . (221) 
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This recursive equation in Fourier space is 



5(2) _ m(m - 2) ~(2) _ m(m + 2) ~(2) 

2 ^m 



167r'^ 



OpnsVi ^T;*^^^ ^^ ^-V^' ^^ ■ ^-V^' — - A o 

ZLUbn n,J'„, , ,2 -^m-1 f™ I 1^2 -^m+1 ~ o "m,2 ; 
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for which we find the solution 



y. 



(2) 



vn 



3 m^ 



■2™-l«'l I cosh h msinh -^ — 

V 2 2 



vahd for m > 2, otherwise 3^m vanishes. 

Next, we consider the equation for Xm- Recall that we are expanding it as A'„ 



(222) 



(223) 



+ 



,2v(i) 



,4 v(2) 



^(1) 



g'^Xm + g^Xm + . . . . Here, Xm is the solution to (209), for which we have only quoted 



^(1) 



(2) 



X^ . This however does not contribute to the equation for Xm ■ For the latter we get the 
usual recurrence equation. In the rhs there are contributions from the solution 3^m , Cm and 
the fermionic convolution. In Fourier space the recurrence equations are 

2 






K2) 



mim - 2) ~(2) _ m(m + 2) ~(2) 



(m + 1) 



2 '^''m+1 



cosh ^3^(2) 



2 ^m 



(224) 



TT 



47r 



rTT (5m,2 



where 3^m is the the solution (223). This equation was also solved using the resolvent 
presented in F.3. We write only X2 



(2) 



X. 



(2) 



1287r 



3 



Li2 (-e-l"'l/2) sinh ^ + Li3 (-e-l'^l/^) cosh f + 1 

167r3 



(225) 



(5e-3l-l/2 _ 9e-l-l/2) 



[2coshf log(e-l"'l/2 + i) _i] . 



27 ' '3 

which is the only component that enters in the fermionic TEA equations. Then, from the 
fermionic Y-functions, we get 



$(2) ^ ^(2) ^ 27r5 



57r 



;-i)"i^n-i + 



32( 



n^ 



^vr(-l)" 



TT 



ir. 



1 



^.^A'i^^-^.^J^f 



(226) 



We compute the convolutions in Fourier space. 



$(2) ^ ^(2) _ _T^3^-H + 87^3 iog(i + e-l^l/2) + 327rLi3(-e- 



hl/2^1 



(227) 



and then 

$(2) + ^(2) 



Att' 



3(1 + m2 



2^' ,1 
+ ^ ( 1 



2tiu 



sinh(27rM) 
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+ 



u* 



1 - 



2tiu 



2tx'u 



2„,2 



sinh(27rn) 



(228) 



Evaluating for u = 2gu we obtain for the leading weak coupling limit 

^^^\2gu) + ¥^\2gu) = ^ + ^ + 0(/) 
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Now we have to compute up to the 0{g'^) convolutions in (94). If we use^^ 
we can re-write the fermionic convolutions in the TBA equation (94) as 



,(10)- 



?(io): 



^ 



2K\"'* log ^ - 26',,"'* log $ = KJ log - + Ka,y* log ^$ 



(229) 



(230) 
(231) 



(232) 



It is important to recall that we need only the m — )■ limit of this. For the order g'^ of the 
first term in (232) we get 

167rV r^ ,n^ 327r V 



-Ka{0) 



3a 



(233) 



The second term in (232) is 

i^,,,ilog^$ = [dvK,^y{0,v)[\og{-A) + 2g^¥'\v) + ^'-'\v)) 



(234) 



-2s 



{¥^\v) + <l>«(t;))' (?^ + 2 (*(2)(^) ^ $(2)(^;)) / + 



where its g"^ order is 



2^4 






Thus, the total (7^ order of (232) is 



32n^g^ SvrV 127rV 



3a 3 

For the remaining convolution in (94), 2 
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-7?(10) , ;o(10) 



* log(l + 3^fe) we use 
<'\Q^ ^) + ^6-2(0, v) = K,.,{v) + Y. ^^^^^ + ^(^') • 



g — 3 
2 



(236) 



(237) 



^Note the difference with Ky^a defined in (200) 



57 



We go to Fourier space, where the term order g^ is 

0-3 

2 

Kb-i + 2_^ Kb^2j 



J 2 



W 



Svr^ _i,„i 167r^a — 1 a> 

p. II P 2 



62 3 3 a 

We Fourier transform back and evaluate for m — )■ and we get 



(238) 



g — 3 

2 



Kb-x^ J2 K, 



b+2j 



y, 



(2) 
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327r2 327r2 Sn"^ 
3a2 3a 3 



(239) 



u=0 



Thus, for Aconv up the 3-loop order we have 



2„2 



A. 



27c^g 



+ 



6tt^ 16tt^ 



45 3a2 



g' + 0{g' 



(240) 



Which, together with (205), gives rise to 

4 



C, = 4(-l)V + 8(-l)'^ 



TT^- 



g'^ + n-ir 



TT^ 47r2 20 



/ + 0(/) . (241) 



F.3 Recurrence resolvent 



For solving perturbatively the small TBA equations we are faced with certain recurrence 
equations, all of same form but with different inhomogeneities. They can be solved in terms 
of the resolvent Xn™ — /(^^) n),for n>m and x^ = f{n, m) for n < m with 



r- (cosh f + n sinh ^ ) (cosh ^ sinh ^^^ — m sinh ^ cosh ^ 

/ \iiLj iLj — '-'V '-' 



mn{coshw — 1)^/^ 



n\w\ 

-e-^- (242) 



This resolvent Xn™ is the solution to 



in T) 

2 cosh -™ -™ 



2^2 



1 An An— 1 An+1 "n 



(243) 



This can be transformed into our recurrence relations by setting Xr. 



^X- 



m^ — l 
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